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PRESENTAZIONE DELLA COLLANA (Presentation of the series) 


Da novembre 2013, con il XXIX ciclo, è attivo presso il Dipartimento di Ingegneria dell’Informa- 
zione dell’Ateneo fiorentino il Dottorato di Ingegneria dell’Informazione (DINFO). Tale dottorato, che 
fa parte della Scuola di Dottorato di Ingegneria dello stesso Ateneo, si propone di formare studiosi in 
grado di svolgere attivitá di ricerca di alto livello in svariati settori dell’Ingegneria dell’Informazione 
e affini, quali Automatica, Ricerca Operativa, Ingegneria Elettrica, Elettronica, Elettromagnetismo, 
Informatica, Bioingegneria e Telecomunicazioni. La presente collana si propone di raccogliere alcu- 
ne monografie legate a corsi specifici del Dottorato, cercando anche di fornire una panoramica delle 
attività di ricerca che si svolgono nell’ambito del Dottorato stesso. Questo volume, in particolare, 
rappresenta un’introduzione alle attività svolte presso il Laboratorio RF, Microonde ed Elettromagne- 
tismo (DINFO) nell’ambito della localizzazione indoor basata su segnali radio, e funge da materiale di 
supporto per il corrispondente corso di dottorato intitolato Introduction to Indoor Wireless Positioning 
Systems, tenuto da Stefano Maddio. 


In November 2013, at the Department of Information Engineering of the University of Florence, 
the Doctorate of Information Engineering has been activated. PhD courses, belonging to the Enginee- 
ring PhD School of the University of Florence, are aimed to the education and training of graduates 
able to perform high level research activities in different areas of Information Engineering and similar, 
as Automation, Operative Research, Electric Engineering, Electronics, Electromagnetics, Information 
Technology, Bioengineering and Telecommunications. This series is aimed to the collection of some 
monographs related to specific PhD courses, trying to provide an overview on the research activities 
carried in the context of the doctorate. This volume, in particular, is an introduction to the activities 
carried out at the Laboratory RF, Microonde ed Elettromagnetismo in the field of indoor wireless loca- 
lization, and serves as a support material for the corresponding doctoral course entitled Introduction 
to Indoor Wireless Positioning Systems, held by Stefano Maddio. 


Firenze, gennaio 2017 


Luigi Chisci 


PREFAZIONE 


Di recente, la capacitä di localizzare persone e beni é diventata la forza trainante di un nuovo 
contesto di applicazioni spazialmente consapevoli (context-aware), che vanno dai servizi di sicurezza e 
protezione, alla navigazione di ambienti per fini commerciali, dalla gestione delle merci nei magazzini 
di stoccaggio, all’intrattenimento. 

Sul fronte degli ambienti aperti, ovvero in un contesto indicato comunemente in letteratura scien- 
tifica come outdoor, lo standard tecnologico ormai consolidato è il sistema Global Positioning System 
(GPS), disponibile oramai su tutti i dispositivi personali. Il sistema GPS richiede per un corretto 
funzionamento la comunicazione con almeno quattro satelliti in orbita, ed è necessario che questa 
comunicazione avvenga in condizioni ottimali di Line of Sight (LoS). La costellazione dei satelliti si 
trova in orbita geostazionaria a circa 20000 km: data la distanza percorsa dai segnali, si comprende 
bene come piccoli errori percentuali sulla stima della distanza siano particolarmente gravosi per la sti- 
ma posizionale dell’utente finale. Per rendersi conto dall’ordine di grandezza, un errore relativo dello 
0.001% sulla stima della distanza utente-satellite corrisponde ad un errore assoluto sulla distanza pari 
a 2 km. 


Figura 1: Caratterizzazione di un'antenna intelligente, cuore del sistema di localizzazione sviluppato pres- 
so il Laboratorio RF, Microonde ed Elettromagnetismo del Dipartimento di Ingegneria dell'Informazione 
dell'Università di Firenze. 


Sebbene siano stati apportati molti miglioramenti nel campo della navigazione satellitare (su tutti si 
ricorda il nuovo sistema Galileo sviluppato dalla comunità europea), l'accuratezza di questa tecnologiaè 
tuttavia fortemente limitata in ambienti chiusi (indoor), dove la ricezione dei segnali dai satelliti è 
deteriorata a causa dell’attenuazione del segnale e degli effetti delle riflessioni multiple, ovvero dal 
multipath. L'ambiente indoor, specie se complesso come un museo o un grande spazio commerciale, è 
uno scenario ostile al sistema satellitare, a causa della difficoltà incontrata dal segnale a propagarsi 
attraverso le opere murarie senza subire distorsioni inaccettabili. Lo scenario indoor non è l’unico 
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caso in cui il la navigazione satellitare è non disponibile o estremamente inaffidabile: ad esempio aree 
all’aperto caratterizzate da un denso fogliame e aree aperte ma estremamente urbanizzate, come i 
moderni centri cittadini affollati di grattacieli, rendono comunque ardua la stima posizionale satellitare. 

Negli ultimi decenni, molti sforzi sono stati dedicati alla sviluppo di sistemi di localizzazione di tipo 
locale, pensati per ambienti indoor e non dipendenti dall’accesso alle costellazioni di satelliti. Questo 
sviluppo indipendente è stato in parte trainato dallo straordinario sviluppo delle reti di sensori — 
Wireless Sensor Network (WSN) - nate per monitorare in maniera autonoma aree più o meno vaste, 
senza l'intervento dell’uomo. Queste reti hanno trovato applicazioni in molti campi interdisciplinari 
e sono spesso interfacciate con le infrastrutture preesistenti, come le Wireless Local Area Network 
(WLAN). In tutti questi campi, l'associazione di una grandezza monitorata ad una precisa posizione 
geografica è pressoché obbligatoria, pena l'inutilitá del dato misurato, che resta privo di contesto. 
Questo approccio si sta evolvendo nell’ancor più dirompente campo dell’ Internet of Things (IoT), e 
di tutte le sue ricadute quali la domotica e l’automazione avanzata. 

È chiaro che, per ovvi motivi economici, un sistema di localizzazione indoor con fini domotici, 
commerciali o accademici, non può competere con una costellazione di satelliti in orbita in termini 
di costi e prestazioni. Inoltre, le caratteristiche di un sistema di localizzazione indoor dipendono 
molto dalla scala geometrica che caratterizza il singolo problema, preso caso per caso. Se lo scopo del 
sistema di localizzazione è assegnare la posizione ad un visitatore in un museo, una precisione sub- 
metrica può ritenersi adeguata, ma se l’applicazione è tele-chirurgica, allora la scala della precisione 
deve essere almeno micrometrica. È chiaro quindi che la complessità, sia software che hardware, e 
di conseguenza i costi che un sistema di localizzazione deve sostenere, possono essere estremamente 
variabili, quantificabili soltanto una volta esplicato il caso concreto. 


Figura 2: Costellazione di sensori appesi al soffitto nei locali dell'Università di Firenze. Questa 
configurazione permette una localizzazione radio non invasiva di nodi mobili e utenti. 


Seguendo questa filosofia, presso il laboratorio RF, Microonde ed Elettromagnetismo del Diparti- 
mento di Ingegneria dell’Informazione dell’Università di Firenze, sono stati sviluppati molti prototipi 
per sistemi di localizzazione indoor basati su segnali radio. L’idea di base è quella di operare con 
un sensore che sia in grado di svolgere la funzione di ancora sia per sistemi di trilaterazione (basati 
sull’intersezione di cerchi), che di triangolazione (basati sull’intersezione di rette). La chiave di questo 
sistema è un nodo wireless, illustrato in Figura 1, dotato di una antenna intelligente in grado di ricon- 
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figurare il proprio fascio. Sulla base di questo nodo elementare, sono state sviluppate reti di sensori 
scalabili, in grado di localizzare un nodo mobile sulla base di comunicazioni radio standard, operanti 
secondo lo standard WiFi 802.11. Un esempio di rete di sensori è mostrata in Figura 2. 

Il presente testo Introduzione ai Sistemi di Localizzazione Indoor vuole essere un primo riferimento 
per chi intende cimentarsi in questo promettente campo di studi in rapida via di espansione. Questo 
volume è infatti il secondo numero della collana dal titolo <Quaderni del Dottorato di Ricerca in 
Ingegneria dell’Informazione dell’Università di Firenze> a supporto dell’omonimo corso di Dottorato 
dell’Università di Firenze. 

Nel caso di interesse di questo testo, le misure coinvolte nella stima posizionale sono esclusivamente 
i segnali a radio frequenza, operanti nella comunicazione tra l’oggetto da localizzare e il sistema stesso, 
in piena compatibilità con il paradigma portato al successo dall’approccio IoT. Nei 3 capitoli di cui si 
compone il presente testo vengono affrontati in modo generale, ma comunque rigoroso, le strategie di 
localizzazione base che fungono da fondamento per i moderni e complessi sistemi di posizionamento 
che si stanno affacciando sul mercato. Viene inoltre fornita una dettagliata bibliografia per permettere 
al lettore l’approfondimento di specifiche tematiche legate a questo mondo. 


Firenze, luglio 2016 


Stefano Maddio 
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1. MODELLI PER LA LOCALIZZAZIONE RADIO 


L’obiettivo di un sistema di localizzazione (o di posizionamento, come verrä indicato alternativa- 
mente in questo testo) & di attribuire ad un mezzo mobile indipendente, indicato come nodo target o 
semplicemente target nel prosieguo di questo testo, una posizione definita rispetto ad un riferimento 
noto, ovvero un insieme di 3 coordinate spaziali che rappresentino in maniera non ambigua la col- 
locazione di quest’ultimo nel riferimento del sistema. Le coordinate spaziali possono essere di tipo 
cartesiano, sferico o cilindrico, avendo opportunamente definito il tipo di riferimento a cui queste sa- 
ranno riferite. E implicito che assieme alla posizione sia indicata una precisione ad essa associata, nel 
senso che, come ogni apparato di misura reale, il sistema di localizzazione è affetto da un’incertezza 
legata alla propagazione dell’errore. 

Un sistema di localizzazione si caratterizza con una precisione adeguata alla scala che il sistema 
stesso si propone: nel caso di un geo-posizionamento satellitare un errore di qualche metro può essere 
adeguato allo scopo; se invece l’idea è quella di individuare la posizione di un visitatore in un museo 
sarà probabilmente indispensabile ricercare un errore sub-metrico. La sequenza storica ordinata delle 
posizioni successive di un target si configura poi come tracciamento (tracking) dell’oggetto, e usufruisce 
di molti altri risultati mutuati da altre discipline, come ad esempio il filtro di Kalman. 

In generale la localizzazione richiede l’osservazione di un sottoinsieme di grandezze fisiche. Questi 
valori di misurazione caratterizzano una certa posizione all’interno dell'ambiente monitorato, pertanto 
è evidente che le osservazioni per scopi di posizionamento devono mostrare una dipendenza spazia- 
le/locale. Questo capitolo è dedicato ai principi di posizionamento, che sono tipicamente applicati 
quando ciò che si osserva e misura è una forma d’onda a radiofrequenze o alle microonde. Benché 
esistano altre grandezze monitorabili per determinare la posizione, quelli basati sullo stesso apparato 
utilizzato per la comunicazione radio hanno l’indubbio vantaggio di non dover duplicare la sensoristica, 
e sono perciò particolarmente adatti al moderno paradigma dell’ Internet ot Things e della domotica. 

Nell’eseguire il compito di localizzazione, con la collaborazione del nodo localizzando, si possono 
individuare quattro fasi: 


1. osservazione del segnale; 

2. estrazione di parametri del segnale dipendenti dalla posizione del target; 
3. elaborazione dei parametri; 

4. post-elaborazione basata eventualmente sullo storico delle osservazioni. 


Un sistema di localizzazione può essere definito effettivo se riesce a individuare la posizione con 
adeguata precisione — per esempio persona in un museo con un errore sub-metrico — ed efficace se 
riesce nel suo intento con risorse minimali (in termini di hardware, software, tempi, ecc). Quest’ultima 
osservazione può decretare il successo o il fallimento di un sistema di localizzazione in un mercato 
estremamente competitivo. 

Dal punto di vista matematico, il modello di osservazione può essere espresso in modo molto 
semplice e sintetico nella seguente espressione: 


r=f(z)+n (1.1) 
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dove r e il vettore delle grandezze misurate, ovvero una specifica caratteristica legata alla grandezza 
fisica utilizzata per rivelare la posizione, x la posizione incognita del mezzo mobile, f(x) è una funzione 
vettoriale nota, eventualmente non lineare, che leghi la posizione x alla grandezza r; infine n è un 
rumore additivo supposto a media nulla. Il tipo di grandezza monitorata determina l’approccio per 
l'elaborazione che poi si traduce nella determinazione della posizione. Si può altresì affermare che la 
relazione f(-) associ la posizione x nel sistema di riferimento noto allo specifico segnale N-dimensionale 
osservato r. Quest’ultima associazione è affetta dal rumore n, che anche in un contesto semplificato 
non può prescindere dal rumore termico insito nel sistema elettronico che determina la grandezza r. 
Il sistema di localizzazione si intende composto da un numero variabile di ancore, ovvero sensori di 
posizione nota, che monitorino il nodo localizzando attraverso la misurazione della grandezza r. Questi 
nodi ancora si intendo posizionati in una fase di setup e costituiscono l’infrastruttura di riferimento. 
Il loro insieme sarà indicato come costellazione. 

Lo scopo del presente testo è l’analisi dei soli approcci basati su misure a radio-frequenza; si 
aggiunga inoltre che il problema verrà affrontato con un approccio bidimensionale, ovvero la posizione 
del nodo localizzando si intenderà come completamente nota una volta note due delle sue 3 componenti 
spaziali, ad esempio x e y nel caso il sistema di riferimento sia cartesiano (come nella maggioranza 
dei casi di interesse). Questa semplificazione è legittimata dal tipo di dominio che un sistema di 
localizzazione radio, come quelli trattati in questo testo, è in grado di coprire. Infatti, non soltanto 
la posizione bidimensionale è compatibile con tutti i casi di interesse pratico che la nuova frontiera 
della localizzazione indoor presenta (si pensi ad esempio al caso di monitoraggio per la fruizione dei 
beni museali, alla navigazione in spazi commerciali, ecc.), ma oltretutto risulta relativamente semplice 
estendere la maggior parte dei risultati al caso tridimensionale effettivo. 

In questo contesto le strategie principali per la localizzazione possono essere divise in 2 macrogrup- 
pi: triangolazione (triangulation) e trilaterazione (trilateration), che a loro volta possono specializzarsi 
in una forma differenziale o in altre implementazioni ancor più specifiche, queste ultime nate in genere 
da esigenze pratiche relative alle grandezze misurate. 


1.1 Trilaterazione 


La trilaterazione, o multilaterazione come viene più propriamente definita nel generico caso in cui 
siano disponibili L riferimenti, stima la posizione del target basandosi sulle proprietà geometriche dei 
triangoli. 

In particolare, la posizione di un’unità mobile si ottiene come risultato di un sistema non lineare 
operante sull’insieme di distanze tra il nodo target e ognuna delle ancore di riferimento. Definita infatti 
la distanza dı dalla prima ancora, la posizione del target è vincolata a giacere su di un cerchio di raggio 
dı dal centro di questa. Data poi una seconda stima da rispetto ad un seconda ancora — e quindi da 
un secondo centro — viene a determinarsi un secondo cerchio su cui il target è vincolato a giacere. 

Due cerchi possono intersecarsi soltanto in due punti, che corrispondono a due possibile posizioni 
del target. Si esprime quindi la posizione del target come r = (x,y) e si suppone nota la posizione di 
due ancore, rispettivamente pari a (11, yi) e (12, y2). La distanza tra le ancore e il target è quindi pari 
a di, d2, e la relazione tra le grandezze in esame risulta: 


(e- 21) + (y — y) = dí 


(z = x2)? + (y — yo)? = di (1.2) 


Come si evince dalla Figura 1.1, la posizione del target è quindi individuata dall’intersezione dei due 
cerchi, ma presenta un’ambiguitä. Per risolverla è necessario ricorrere a una terza misura rispetto ad 
una terza ancora di posizione (x3, y3), che individua una distanza dz con il target. In questo modo il 
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Figura 1.1: ambiguita della posizione del target con due misure di distanza. 


sistema si traduce in: 
(2-21)? +(y-y) =d 
(x — 12)? + (y — ya)? = di (1.3) 
(x — 23)? + (y — ya)” = di 


che ammette quindi una soluzione univoca. 

L’approccio classico per la soluzione di questo sistema si basa sul teorema della corda. A tal scopo si 
supponga, senza perdere di generalitá, di poter fissare il centro del riferimento assoluto nella posizione 
dell’ancora #1, e che l’asse x sia quello che la congiunge all'ancora #2. Questo implica che x] = 0, 
yı =0e z2 = D, y = 0. Con riferimento alla Figura 1.1 si ha che 


x + y? = dî iù 
(€ — DP + y? = d da 
y = ag 
sottraendo la seconda equazione alla prima si ottiene 
x? — (1 — DP =22D- D’? = dí — dí (1.5) 
da cui, facilmente 
aD, d-d 
a 2: T 3D 2 (1 6) 
T=+Vd;° — ?? 


dove il segno + è dovuto all’ambiguità precedentemente descritta. 

Per risolvere questa ambiguità è sufficiente aggiungere una terza ancora, come indicato in Figu- 
ra 1.2. In questo caso infatti ogni coppia di cerchi produrrà due stime posizionali formalmente identiche 
alle (1.6), ovvero sei stime in tutto. Di queste, almeno tre saranno coincidenti, e corrisponderanno alla 
reale posizione (x, y) del target T. Si può pertanto concludere che, nel caso ideale, il numero minimo di 
ancore necessarie per trilaterare la posizione è pari a 3. Tuttavia, a causa di inevitabili errori di misura 
nella stima delle distanze di, da e da il sistema di trilaterazione a tre ancore è risolvibile esattamente 
soltanto in via teorica, perché l’effetto del rumore è di perturbare i tre cerchi, impedendo l’esatta in- 
tersezione e quindi la determinazione della posizione del target. Si può affermare che anziché con una 


INTRODUZIONE AI SISTEMI DI LOCALIZZAZIONE INDOOR 


i 


A 
x 
>, 
, 
. 
. 
. 
. 
. 
s 
— ia 
> 
; 
; 
2 
; 
È 
© 
s 
. 
, 
, 
, 
, 
y 
A 
N 
i 
i 
N 
E 
ili 
} 
i 
P 
; 
; 
; 
; 
; 
; 
y 
; 


ra 
a 
< 
WY 
en. 
~ 
s 
< 
a 


` 


a 


Figura 1.2: Posizione individuate in modo non ambiguo mediante un sistema di tre ancore 


distanza deterministica si ha a che fare con un intervallo di possibili distanze, dal momento che nel 
caso pratico € almeno possibile stimare un intervallo di fiducia sulla misura di d, individuando quindi 
gli estremi entro cui d si trova con sicurezza statistica. Si parla talvolta di pseudo-range, specialmente 
nelle implementazioni basate sul calcolo dei tempi come misura indiretta per la distanza. 
Dall’ispezione diretta si può inoltre comprendere come l’errore con cui è affetta la stima delle 
distanze d possa avere effetti diversi in conseguenza della posizione localizzanda. L’accuratezza del- 
l’espressione an che determina la coordinata x del target, dipende molto dal rapporto tra d e D. 
Se quest’ultimo e molo piccolo, l’errore su d verrà esaltato nel calcolo causando una propagazione 
nefasta dell’errore. In Figura 1.3 questo concetto è mostrato in modo grafico. A sinistra, l’intersezione 
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Figura 1.3: Propagazione dell'errore per effetto della posizione delle ancore rispetto al target 


di due cerchi di uguale raggio si modifica in seguito ad una perturbazione sulla misura degli stessi (per 
semplicità l’errore è supposto identico su entrambi). In seguito a questa deformazione, la posizione 
del localizzando è incerta all’interno dell’area sovrapposta dei due segmenti circolari. In altre parole 
la posizione del target potrebbe essere ovunque all’interno della regione grigia 

A destra, la medesima situazione si ripete nel caso in cui la distanza tra le ancore sia molto 
minore di quella nel caso a sinistra. A fronte di una medesima perturbazione sui raggi dei cerchi, 
l’area identificata dall’intersezione dei segmenti circolari è maggiore che nel primo caso e molto più 
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deformata in un senso rispetto all'altro. Questa situazione avrá quindi un peso maggiore sull'incertezza 
della stima finale. 

E abbastanza chiaro che, a prescindere dalla perturbazione dovuta all’errore, la stima posizionale 
migliora se l’intersezione viene eseguita a partire da piú coppie di cerchi e quindi pit distanze, poiché 
questa sovradeterminazione raffina l'intersezione tra gli intervalli di sicurezza con cui le grandezze d 
sono stimate. Tuttavia, anche in presenza di un sistema sovradimensionato, la soluzione della trila- 
terazione in maniera diretta come intersezione di cerchi non è generalmente il modo migliore con cui 
procedere. La migliore strategia passa infatti attraverso la riformulazione del problema in maniera pit 
strutturata, in modo da sfruttare anche i risultati provenienti da altri ambiti di ricerca. 

Nel caso generale di una costellazione di L ancore, il modello di misura in (1.3) si riformula in 
maniera estesa nel seguente modo: 


d= ||e - aıl| = V(e- z1)? +(y-y), 1=1,2,...,L. (1.7) 


Pertanto, la !-ima misura della distanza indicata con rr, si può indicare con 


TRI = di +nri = V(x = x1)? + (y vi)? + NR; 1 RP Y (1.8) 


che può essere a sua volta espressa in forma vettoriale come: 


rri=fr(x)+nR (1.9) 

dove 
TR =(re,vraa rr), (1.10) 
NR, = (nea NR) nr], (1.11) 


La funzione vettoriale fr(x) è quindi composta da componenti parametrizzate in x, pertanto il 
problema della localizzazione si traduce nella stima di x data la grandezza misurata rp, che si intende 
sempre affetta da errore come espresso dalla relazione in (1.11). 

Riscrivendo il tutto esplicitamente, la relazione in questione diventa: 


TRA Y (a - 21)? += y)? NR,1 
RAS y (1 — 22) n (y — y2) m Bo (112) 
rail Ivan lai 


Un modello alternativo per questo tipo di osservazione è quello basato sull’approccio probabilistico. 
In questo senso, la generica osservazione eseguita dall’ancora /-sima rispetto al target è vista come 
una densità di probabilità di tipo normale (gaussiana), esprimibile come: 


1 1 
p(rrRı) = exp (rr a) (1.13) 
A ¡2707 | | 20k 1 


La media di questa distribuzione è ovviamente d; e la deviazione standard dovuta all’errore è pari a 
or. In questa accezione, il modello (1.12) può essere riscritto in forma di statistica L-dimensionale 
nel seguente modo: 

1 


e 
(27)"’? |Cp]1/2 


plen) = xp (¿rn a)" Ctra ~ a) (1.14) 
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dove la matrice di correlazione Cr può essere espressa come 
Cr=E[(rr-d)(rr- d)'], (1.15) 


essendo l’operatore E [a] il valore atteso della grandezza a. Nel caso in cui le distribuzioni di probabilità 
siano indipendenti, ovvero nel caso di rumori di misura scorrelati tra le varie ancore, la matrice di 
correlazione può essere semplificata: 


Cpr = E [ngng] = diag (0% ,,0% 2) ---,0h,L) (1.16) 


E di conseguenza anche il modello (1.13) assume la forma di una distribuzione gaussiana multidimen- 
sionale: 


p(rr) ~ N (d, diag(0r,1,0R,2,..., 9R,L)) (1.17) 
che puo essere espressa nel seguente modo: 
1 1 (rr — di)? 
p(rr) = exp (1.18) 
(20) P TH ory 2 > el 


1.1.1 Stima della distanza basata sul tempo di volo 


Gli approcci basati sulla trilaterazione sono stati tra i primi ad essere impiegati per la stima 
della posizione del target incognito, sfruttando misurazioni temporali dei segnali radio. Queste ultime 
possono infatti essere facilmente convertite in misure di distanza conoscendo la velocità dei segnali in 
gioco, [1-3]. 

Il tempo impiegato dal segnale per percorrere lo spazio che separa il target localizzando posizionato 
in r dalle singole ancore è tipicamente definito tempo di volo, o con una diffusa nomenclatura inglese 
Time of Flight (ToF); nome alternativo è quello di tempo di arrivo, ovvero Time of Arrival (TOA) 
in Inglese, che enfatizza invece l’istante in cui il segnale viene rilevato. La velocità del segnale è 
evidentemente la velocità della luce nel mezzo di propagazione, tipicamente l’aria. 

In genere si assume che ogni nodo ancora abbia accesso allo stesso clock condiviso dall’intera 
costellazione. Si indica dunque con t; il tempo di arrivo alla i-esima ancora. Assumendo che il segnale 
sia trasmesso al tempo tm, sia 7; il ToF tra il target e la generica ancora (x;, yi), distante d; dal target. 
Utilizzando queste informazioni, è possibile scrivere che 


d; i 2 i 2 
e WE (1.19) 
Cc C 


che assume il ruolo della (1.7). Senza perdere di generalitá, si pud supporre che la sorgente — il target 
- emetta il segnale al tempo zero, e che il sensore nell'ancora I!" lo riceva al tempo tı; in questa 
situazione il modello di misura in (1.3) si riformula in maniera estesa nel seguente: 


d 1 


TTOA 1 = č = „vie xi)? | (y y? | NToA,l; l= 152,232; (1.20) 


Di conseguenza i modelli per trattare questo caso sono identici a quelli in (1.9) e (1.13) a parte il 
fattore di scala dovuto alla velocitá della luce. 

Quest'ultima è appunto la principale difficoltà che affligge i sistemi basati sul principio ToA, dato 
che l’elevatissimo valore di 3.00 x 10° m/s, rapportato alle distanze metriche degli ambienti indoor, 
rende necessarie precisioni e coerenza estreme. Per ottenere con un valore di precisione adeguato il 
ritardo di propagazione del segnale, i riferimenti temporali di tutte le ancore e del target devono 
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essere uguali, ma questo requisito è difficile da raggiungere. Dal lato della costellazione di ancore è 
possibile impiegare una ragionevole quantità di risorse per sincronizzare le basi tempi, utilizzando 
orologi altamente stabili. L'eventuale asincronia tra la ancore può inoltre essere misurata e fornire così 
il dato occorrente per per compensare tale deviazione. Tuttavia, il problema di differenti basi temporali 
tra la costellazione di ancore ed il nodo target permane, almeno per dispositivi commerciali di massa. 
Questi ultimi sono infatti costruiti con componenti a basso costo e possono al massimo fornire una 
ragionevole stabilità a breve termine. Questa incertezza implica un’ineliminabile perturbazione sui 
cerchi, causando un’intersezione rumorosa come quella indicata in Figura 1.3. 


1.1.2 Stima della distanza basata sulla potenza — RSSI 


Con il termine Received Signal Strength Indicator, o RSSI, si intende un valore numerico che indica 
la potenza ricevuta dai sensori di riferimento durante una comunicazione con il target, [4, 5]. Si assume 
comunemente che la potenza ricevuta segua il modello di decadimento polinomiale, tipico del modello 
di Friis, [6, 7]. Tale potenza è infatti supposta dipendente quella trasmessa, dalle costanti di perdita 
lungo il percorso — path loss constant — e dalla distanza tra sorgente e sensore. Il fatto che ogni 
trasmettitore conforme agli standard IEEE sia nativamente in grado di misurare il valore dell’RSSI 
dei pacchetti radio ha reso molto popolari i metodi basati su questo approccio, [8-13]. 

Assunta una potenza trasmessa P;, in assenza di disturbi la potenza media ricevuta dall’° sensore, 
indicata come P, ¡, è modellata da [14]: 


Ppa = Ki Pied, * = Kı Pæ — æl, 1,2,...,L (1.21) 


dove i coefficienti K, tengono conto di tutti quei fattori che influenzano la potenza ricevuta ma che 
non dipendono dalla distanza, come ad esempio il guadagno delle antenne usate per la comunicazione 
(che pertanto dovranno essere ben caratterizzate). Il coefficiente a è invece la costante di path-loss, 
che in funzione dello scenario può essere un numero variabile tra 2, valido per lo spazio libero, a 5 o 
6 in un ambiente chiuso, a causa di riflessioni e altri fenomeni diffusivi di ordine superiore. Esistono 
molti modelli empirici ed euristici che permettono di conservare la validità del modello a patto di 
aumentare appunto il valore di tale costante. Inoltre, svariati anni di sperimentazioni hanno messo in 
evidenza come il rumore additivo alla grandezza RSSI sia di tipo lognormale — dove con tale termine 
si intende una distribuzione di probabilità di una variabile aleatoria r il cui logaritmo In (r) segue una 
distribuzione normale, [15]. 

Basandosi quindi sulle precedenti considerazioni, il modello di osservazione completo può essere in 
questo caso espresso come: 


In Pri =lnK1+lnP,- alnd + NRSSIl; E RA Y (1.22) 


dove la distribuzione nrsgr,ı è normale. Per semplicità si assume che nrss1,1, | = 1,2,...,L siano 
processi gaussiani a media nulla caratterizzati da varianze {c}gg,}. Definito il valore 


rrssıı=mP.ı-nK, — ln P, (1.23) 
il modello si semplifica in: 
rrssıı= —alndı +nrssı., 1=1,2,...,L (1.24) 
che assume la seguente forma vettoriale: 


rrssi = frssi(x) + nessi (1.25) 
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dove 


rrssi = [Rssı.1,TRsst2,..-TRssı,r) (1.26) 


NRSSIT Messia; NRSSI,2) +++ nessi.) ; (1.27) 


In y (z — 21)? + (y = ya)? 


In y (£ — 292)? + (y — ya)? 
Frssile) = p = ~a i i ee) (1.28) 


n/a 0? 


Pertanto, il problema della localizzazione del target basato su stime della potenza ricevuta può essere 
formalizzato come la stima di x dato il set di misure formalizzate come rrgssr. 

A livello pratico lo schema di posizionamento è più semplice di quello usato per le misure dei 
tempi di volo, dal momento che non sono necessarie sincronizzazioni. Tuttavia, la complessità del 
canale radio implica delle perturbazioni notevoli (e spesso con distribuzione non uniforme) sulla stima 
della potenza ricevuta. 

In questa implementazione, la formulazione basata sulla densità di probabilità assume la forma: 


1 1 
p(rrssI) = exp ( (rassi — d) Caggr(ra — a)) = 
(2n)"/? [E p|1/2 2 RSSI 
L 
1 1 Indi)? 
va exp (-3), asst end | (199) 
(27) “TIE \0RSSII os ORSST, 
dove la matrice di correlazione Crggy risulta pari a diag (orssı 1, ORss12,---,ORssı,L) (?)- 


1.2 Trilaterazione differenziale 


La trilaterazione differenziale si basa sulla differenza tra le distanze, ovvero sulla combinazione tra 
i range stimati dalle varie ancore, avendone assunto uno come riferimento, [16]. Con riferimento alla 
stessa distribuzione di ancore di Figura 1.1, la distanza d fa stavolta da riferimento, per cui: 


=. 2)? + (y- y2)? - /(@— a1)? + (yy)? 
dz — dı = y (a — 23)? + (y — ya)? - y (a - 21)? +(y-yı)? 


Geometricamente, la prima relazione definisce un'iperbole con i fuochi in (x1,y1) e (2,42), come 
si puo osservare in Figura 1.4; la posizione del target e ottenuta dall'intersezione di almeno due 
iperboli, quindi in una costellazione che dispone di almeno tre ancore, come illustrato in Figura 1.5. 
Nel caso generale di più ancore si avranno ovviamente più iperboli che a causa del rumore di misura 
non necessariamente si intersecheranno in un singolo punto, andando così a rendere la localizzazione 
ambigua. 

Date L ancore, esistono L(L — 1)/2 differenze di range distinte considerando ogni possibile coppia, 
ovvero: 


(1.30) 


dkı =dk— dı, k,l=1,2,...,Lconk>l. (1.31) 


Tuttavia, soltanto (L— 1) di queste differenze sono indipendenti. Al fine di ridurre la complessità senza 
sacrificare le prestazioni, è preferibile misurare tutte le L(L—1)/2 per poi convertirle nel sottoinsieme 
delle L — 1 stime non ridondanti [17]. 
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Figura 1.4: Luogo dei punti in cui può trovarsi il target quando è nota soltanto la differenza də — dı. 
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Figura 1.5: Posizione individuate grazie all'intersezione di due iperboli ottenute mediante tre ancore. 


Senza addentrarsi troppo nei dettagli di questa strategia, e senza perdere di generalitá, si assume 
che il sensore 1 faccia da riferimento per l'implementazione diretta, per cui la singola osservazione su 


cui si fonda il modello risulta: 


dı =d — dı = y (£ — 11)? + (y — u)? — Y (2—-20)2+(y—y)?, 1=2 


Assunto ciò, la funzione di osservazione vettoriale comprensiva di rumore diventa 


ror = y(x — z)? + (y - y)? — y(2- z1) +(y- y1)? + nori, 1=2 


pee 


L. (1.32) 


„L (1.33) 
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che concisamente in forma vettoriale si scrive: 


TDR = fori(2) + NDR (1.34) 


dove 
TDR = [FDR 2;TDR3;---;"DR 1) (1.35) 


e inoltre 


y (a = 22} + (y — 42)? - y (a — 21)? + (y — 41) 


fonte) = E eat) ui 


Je- +a jean 


Il problema della localizzazione con approccio di tipo triangolazione differenziale si enuncia formal- 
mente come la stima di æ nell’equazione (1.34) dato il vettore rpr. 

Per quanto riguarda l’approccio basato sulla distribuzione di probabilità, in questo caso è necessario 
porre attenzione alla natura dell’errore. Assunto infatti che il vettore di errore npr sia gaussiano a 
media nulla, la pdf associata a questo modello di osservazione si può scrivere: 


1 
(27)? ¡Cor 


(er) = A) (1.37) 


ma in questo caso, la matrice di correlazione Cpr non è diagonale, dal momento che la scelta di 
riferire tutte le stime di distanza alla stima del sensore l = 1 rende i contributi npr, correlati tra 
loro [18]. 


1.2.1 Differenza dei tempi di volo 


Il modo più diretto per implementare la trilaterazione differenziale è attraverso l’applicazione della 
differenza dei tempi di volo, o Time Difference of Arrival (TDoA) sfruttando la dizione con cui è 
nota a livello internazionale. Tale metodo si basa su una misura diretta della differenza dei tempi 
di percorrenza del segnale nel suo tragitto dal target alle ancore. Il grande vantaggio del metodo 
TDoA rispetto al ToA consiste nel rilassamento sulle richieste di sincronizzazione, dal momento che in 
questo caso non è necessario che essa sia assoluta. Si può infatti affermare che il modello differenziale 
si è affermato proprio per permettere di rilassare la richiesta di sincronizzazione assoluta, la cui 
implementazione può risultare estremamente complicata e costosa a livello hardware. 

Sulla basa di questa filosofia è basato il sistema LORAN (LOng RAnge Navigation), utilizzato in 
passato da molte nazioni, quali gli Stati Uniti, il Giappone, l’India e molti paese europei, [19]. Questo 
è un sistema di radionavigazione terrestre tramite onde radio in banda LF, con frequenze da 90 a 
110 KHz - che sfrutta la misura dell’intervallo di tempo tra i segnali ricevuti da tre o più stazioni per 
determinare la posizione di una nave o di un aereo. 

Nell’ambito di questo tipo di sistemi, si assume che la sorgente emetta il segnale al tempo to 
— di per se non noto al sistema — e che il sensore /-simo lo riceva al tempo tı, L = 1,2,...,L con 
L > 3. Utilizzando questo sistema, non è necessaria una sincronizzazione assoluta, dal momento che è 
sufficiente un contatore preciso nel target, ma non obbligatoriamente sincrono con il riferimento delle 
boe. Queste ultime invece debbono essere sincronizzate tra loro, ma questa operazione è senza dubbio 
più semplice rispetto al garantire un target sincrono in maniera assoluta. Tipicamente le ancore sono 
infatti più sofisticate del target, ed inoltre, facendo capo ad un nodo supervisore che è parte della 
stessa rete, sono facilmente sincronizzabili in una fase di setup. 

Il modo più semplice per implementare la strategia basata sulla differenza dei tempi di voli consiste 
infatti nel tenere il target in ascolto, e far trasmettere a tutta la costellazione di ancore il proprio segnale 
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in modo sincrono tra loro. Al target sará sufficiente registrare gli arrivi dei vari segnali dalle ancore 
per costruire una tabella riferita ad un proprio tempo locale. Da tali dati il target è quindi in grado 
di costruire le L — 1 misurazioni necessarie per l'implementazione del metodo. 


1.3 Triangolazione 


La triangolazione & l’approccio duale a quello della trilaterazione. In questo caso la stima della 
posizione del target mobile avviene confrontando la direzione che esso assume rispetto ai riferimenti 
locali delle ancore. Si suppone che queste ultime siano coordinate tra di loro, condividendo quindi lo 
stesso riferimento assoluto. 

La direzione che unisce ancora e target viene spesso indicata come direzione di arrivo, in genere 
con la dizione inglese di Direction of Arrival (DoA), o talvolta anche come Angle of Arrival (AoA). 
Il vantaggio dei metodi basati sulla triangolazione è che non necessitano della stima dei tempi di 
volo (e conseguentemente nemmeno di sincronizzazioni), e non sono neppure basati su stime assolute 
di potenza, come già visto spesso complicate da fenomeni di ordine superiore e dalla necessità di 
calibrazione dei sensori. 

Con dizione mutuata dalla tradizione nautica, si suole indicare come Line of Bearing (LoB) la retta 
che congiunge il nodo target all’ancora di riferimento: di fatto questa nomenclatura è equivalente a 
Line of Sight (LoS), ma quest’ultima viene utilizzata anche in altri ambiti. L’angolo di arrivo è quindi 
quell’angolo & che la retta LoB identifica con l’asse x del target (scelto generalmente come zero di 
riferimento per gli angoli); esso può essere stimato sia dal target rispetto alla costellazione di ancore 
o viceversa da parte delle ancore che osservano il nodo localizzando. In entrambi casi il modello di 
osservazione non cambia, ed è formalmente identificato come segue. 


(x1,y1) (x2,y2) 


Figura 1.6: Posizione del Target T individuata attraverso l'intersezione di due rette tracciate da due ancore. 


Con riferimento alla figura Figura 1.6, siano @1 e da gli angoli di arrivo individuati dal target T 
rispettivamente con l’ancora 1 e 2; la relazione tra le coordinate spaziali dell’ancora e quelle del target 
risulta: 


y-yı = tan (¢1)(x — 21) (1.38a) 


y — y2 = tan (p2)(1 — 22). (1.38b) 


Senza perdere in generalitá, si puó considerare un riferimento locale per cui yı = ya = 0, e sottraendo 
membro a membro le due (1.38) si ottiene la coordinata x: 


_ tandızı — tan daro 


tan dı — tan da ae) 
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Sostituendo questo valore in una delle due (1.38) si perviene alla stima della coordinata y: 


tan dı tan da 
= ; 1.4 
y ee 22) 0) 


Anche in questo caso, come per l’approccio basato sulla trilaterazione, la presenza del rumore di misura 
puo pregiudicare la stima. E infatti possibile che, in presenza di piú di tre ancore, le relative rette 
non si incontrino in un solo punto, ma individuino un insieme di intersezioni. Mutuando il concetto 
di intervallo di distanze introdotto per la trilaterazione reale, in questo caso si ha a che fare con un 
insieme di rette che determinano una sorta di cono entro cui è possibile affermare si trovi l’effettiva 
Line of Bearing. 

Anche in questo caso, inoltre, l’errore di misura affligge la stima di posizioni differenti in modo 
differente. Dalle formule in (1.39) e (1.40) si osserva che in entrambi i casi si ha a denominatore la 
differenza tra due stime angolari, pertanto una stessa variazione può modificare in modo sostanziale 
la stima di x e y a seconda dell’effettivo valore degli angoli. 

Si procede pertanto alla riformulazione del problema in termini più generali, supponendo di poter 
operare con un insieme sovradimensionato di ancore come per il caso della trilaterazione. Sia @; la 
direzione tra il target e l’ancora /-sima, per cui: 


tan oy (222), eat (1.41) 


T — Tl 


In presenza di rumore, la stima dell'angolo formato dalla retta LoB può essere formulato nel seguente 
modo: 


raı=dıtnaı=tan! (£ = n) Frag. l=1,2,...,L (1.42) 
— 21 
ovvero, in forma vettoriale: 
Tan =fa(c) +na (1.43) 
dove 
ra = [ra TA, rA, L] (1.44) 
na = [na nA, naL]; (1.45) 
e infine 
tan Soe 
LI 
tan”! | 2 
LI—XLa 


Salz) = ¢ = (1.46) 


tan”! (222) 
LLL 
Anche in questo caso il problema può essere considerato in modo statistico, con una distribuzione di 
probabilitá simile a quella in (1.13). 


1.3.1 Stima dell’angolo di arrivo basata su array di antenne 


L’angolo di arrivo di un segnale radio coincide con l’angolo di arrivo del fronte d'onda cui esso 
appartiene. Se si assume una propagazione in linea di vista, ovvero ignorando in prima approssimazione 
effetti del secondo ordine come diffrazione, riflessione, etc., € facile comprendere come un array di 
antenne sia adeguato per stimare l’angolo di arrivo nel caso 2D, [20, 21]. 

Con riferimento a Figura 1.7, si assume che un nodo ancora sia equipaggiato con un array lineare 
di antenne — Uniform Linear Array, ULA. La distanza inter-elemento, costante lungo tutto l’array, è 
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Fo r Ip rt 
Figura 1.7: Stime dell’angolo di arrivo basata su un array lineare di antenne. 


indicata con d. Il fronte d’onda in arrivo forma un ¢ con il riferimento dell’array, il che implica che 
ogni antenna riceva un segnale “ritardato” di distanza incrementale A pari a dsin ¢. 

Nell’ipotesi di operare in regime sinusoidale, si puó affermare che i segnali r; ricevuti dagli elementi 
dell’array sperimentano una progressione di fase proporzionale all’angolo di arrivo ¢: 


bp == ( = =) kdsin $ (1.47) 


dove k = = e il numero d’onda dell’array e Ay = Fo é la lunghezza d’onda nel vuoto associata al 
segnale ricevuto alla frequenza fo. 
Ignorando il rumore, il segnale ricevuto dall’elemento /-esimo del set di antenne può essere espresso 


in notazione complessa come [22]: 


ri = spe i(!- 7 )kdsing (1.48) 


dove s; è l’ampiezza del segnale ricevuto. Quest’ultima può essere supposta costante su tutti gli 
elementi dell’array dal momento che la differenza di cammino dsin d è in genere trascurabile rispetto 
alla distanza percorsa dal segnale, e considerato che il modulo del singolo segnale è proporzionale 
all’inverso della distanza. Ovviamente non può essere trascurata nella fase, che invece sperimenta una 
variazione rapida lungo l’array, e che costituisce la chiave per la determinazione dell’angolo di arrivo. 

Con riferimento alla Figura 1.8, i segnali r; vengono pesati da dei fattori a; e quindi sommati. 
Il vettore a = ao, 41,...@7,-1] viene indicato con il nome di steering vector. A seconda della forma 
di questo elemento, il segnale risultante r = va: arri può indicare quale sia la direzione di arrivo 
@ in maniera plurivoca, consentendo quindi moltissime e varie strategie di stima della DoA, [20-22]. 
Sistemi ancor più avanzati sono stati proposti negli ultimi anni, in cui l’array di antenne piuttosto che 
lineare ha una forma circolare. Questi approcci permettono di utilizzare anche le informazioni relative 
al modulo dei segnali incidenti, aprendo la strada a molti altri algoritmi di stima della DoA. 

Si consideri inoltre che operando con queste geometrie è possibile determinare la direzione di arrivo 
anche utilizzando il solo modulo, anziché la sola fase. È questo il caso delle antenne intelligenti, in 
Inglese Smart Antennas, capaci di riconfigurare il proprio fascio in più direzioni, [23-25]. 
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Figura 1.8: Stime dell'angolo di arrivo basata su un array lineare di antenne. 


1.4 Triangolazione differenziale 


Può accadere che il target sia equipaggiato con un sensore in grado di determinare la direzione 
di arrivo del segnale rispetto al proprio riferimento, ma che questo riferimento non sia condiviso 
con la costellazione delle ancore [22]. In queste condizioni non è possibile quindi fissare una direzione 
univoca rispetto alla singola ancora, e conseguentemente non si può procedere alla determinazione della 
posizione tramite l’intersezione delle rette. Stimando l’angolo di arrivo individuato rispetto ad una 


Figura 1.9: Luogo dei punti in cui può giacere il target data la differenza di angoli di arrivo da 2 ancore. 


seconda ancora, il target è tuttavia in grado di fissare in modo certo almeno la differenza degli angoli 
di arrivo rispetto alle due. Questa situazione è duale rispetto a quello che accade nella trilaterazione 
differenziale, in cui è la differenza di distanze a essere nota piuttosto che l’insieme delle distanze stesse. 
È facile dimostrare che, stimato Ad rispetto a due ancore, il target giace senza dubbio sul luogo dei 
punti che osserva il segmento individuato dalle due ancore sotto lo stesso angolo, ovvero un cerchio 
che passa per i due riferimenti in questione. 

Il modo più sicuro di procedere passa dalla riformulazione del problema in termini più generali; 
si supponga quindi di poter operare con un insieme sovradimensionato di ancore, come per il caso 
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. 
.. 


(xiyi) (X2,Y2) 


Figura 1.10: Posiziona del Target T individuata attraverso l'intersezione di due cerchi identificati dalle 
due differenze angolari Adı,2 e Ad» 3. 


della trilaterazione. Si assuma che il sensore l = 1 faccia da riferimento, e si indichi con @; l’angolo 


identificato tra il target e l’ancora l-sima. Il modello di osservazione per il singolo caso risulta di 
conseguenza: 


$= XI 


Abi = di - ¢1 = tant! (+) — tan”! (En) l=2,3,..., L, (1.49) 
XT ZI 


che in presenza di rumore diviene 


roas = Adia + pas = tan! (FE) — tant (HB) pas li=2,3;00aa Le (1.50) 
i LL] 2-2 


ovvero, in forma vettoriale: 
rpa, = fpa(12) +nDA (1.51) 


con il vettore di misure e il vettore di rumore definiti da 


rpa = [rpa.1,"Da2,:- "DA," (1.52) 


NDA = [NDA 1,NDA2)---NDAL]'; (1.53) 


e infine, la funzione che lega misurazioni e posizione: 


—1 / yzy \ _ —1 [ y= 
tan en tan es 
tano! ( 2=%8 \ tan! (u 
T—T3 T—X1 
fpa(a) = Ad = 


(1.54) 


tan! (Eu) _ tan”? (eu) 
LT—LL T—L1 


ritrovando cosi la familiare forma (1.1). 


16 INTRODUZIONE AI SISTEMI DI LOCALIZZAZIONE INDOOR 


Formulazione alternativa 


Esiste infine un modo alternativo per trattare la triangolazione differenziale. Sfruttando il teorema 
di Carnot, altrimenti noto come legge dei coseni, si può affermare che, fissando la prima stima dı 
come DoA di riferimento, vale la seguente uguaglianza per la boa l-sima: 


dî, = dî + dj — 2dıdı cos (1 — (1) (1.55) 


Questa relazione si traduce nel set di equazioni non lineari quando ogni DoA e riferita alla stima 
dell’ancora l = 1: 


(1 — 01)? + (y1 — yi)? = (z — 21)? + (y — y1)? + (2 — 21)? + (y — y)” 
2 (a — x1)? + (y — y1)? y (a — 21)? + (y — yıl)?cos(Adıı) 1=2,3,...,L (1.56) 


da cui si può infine ricavare la relazione: 


(a + (y-n)? + (2—21)? + (y Ha)? = (21 — 21)? = (y1 — w)? 
Adı ı = arccos ; 
1,1 = arccos ( NCEE N a+ Y ) (1.57) 


e quindi: 
(z—z1)?+(y—y1)?+(z—z2)?+(y—y2)?-—(d1,2)? 
arccos 2 (2 02H y)? (e 22)2+(y—y2)? 
arccos | C772 +W=n)+e@=as)’+Hy-ys)’-(dı,3)” 
fpa(x) = Ad = 2 (0214-91). (2-23)? +(y-y1)? (1.58) 


arccos (@—21)? +(y—y1)? +(@— a2)? +(y—yt)?— (dit)? 
2 (221)2+(y=y1)2 Y(@-2L)?+(y-yr)? 


Ciò fatto si può arrivare ad una formulazione completamente polinomiale del problema, evitando quindi 
le funzioni inverse trigonometriche che possono mal-condizionare i procedimenti risolutivi, elevando a 
quadrato i due membri di (1.55) ed ottenendo cosi: 


4 ((& — 21) + (y - yı)*) ((& y + lg yu?) (cos (Adı1))” > 
(@- 1)? +-+ (oa + (y)? (aa? (mm)? 1=2,3,...,L (1.59) 


Date almeno tre stime indipendenti dell’angolo di arrivo da tre diverse ancore, si ottengono due 
equazioni polinomiali indipendenti con il massimo grado di quattro sia in x che in y. Cosi come 
due cerchi si intersecano al massimo in due punti, in questo caso le intersezioni sono 12: è pertanto 
necessario sovradimensionare il sistema in modo che la soluzione reale emerga dall’intersezione delle 
varie curve ridondanti. 
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Gli approcci per risolvere il sistema rappresentato dall’equazione (1.1), e quindi le sue specializza- 
zioni nei vari metodi visti nel capitolo precedente, si suddividono in due grandi categorie: gli approcci 
lineari e i non lineari. 

Per la precisione, quelli lineari sono piú propriamente approcci linearizzati, ottenuti sviluppando 
l’equazione di partenza assumendo alcune ipotesi che dipendono dallo specifico approccio di localizza- 
zione. Ogni metodo linearizzato necessita di attenzioni diverse ed in genere non è univoco il modo in 
cui dall’equazione (1.1) si perviene a un modello lineare, per cui anche l’efficacia dell’approccio dipende 
molto dalla filosofia con cui si opera e dalla perizia con cui si adoperano le informazioni disponibili. 

Per quel che riguarda i metodi non lineari, questi cercano di risolvere l’equazione (1.1) in modo 
diretto, cercando di determinare l’incognita a minimizzando il residuo, costruito a partire dalla spe- 
cifica forma che assume l’equazione (1.1) caso per caso. La funzione di costo viene costruita a partire 
dal residuo implementando un approccio ai minimi quadrati — Least Square (LS) — o la la sua va- 
riante più raffinata con pesaggio — Weighted Least Square, WLS. In generale si indicherà il residuo nel 
seguente modo: 

enr =r-f(®) (2.1) 


dove 2 = [7,7] ' è la variabile di ottimizzazione per æ, che corrisponde allo specifico stimatore. La stima 
del metodo di localizzazione corrisponderà a quel valore di € che rende il residuo il più piccolo possibile 
— da cui la dizione e ad indicare che si tratta di un errore. In questo senso, gli approcci non lineari 
sono più generali di quelli lineari e possono sfruttare i notevoli risultati della teoria dell’ottimizzazione 
apparsi in letteratura scientifica negli ultimi anni. 
Sul fronte delle tecniche lineari, si tratta di convertire l’equazione (1.1) in un set di equazioni 
esprimibili nella seguente forma: 
b= Ax+q (2.2) 


dove b e A sono le grandezze disponibili, derivanti dalla misura, mentre q è il vettore del rumo- 
re di misura trasformato dall’operazione di linearizzazione. Quest'ultimo può essere estremamente 
malcondizionato nella trasformazione, pertanto grande cura deve essere usata in questa classe di 
algoritmi. 

A partire dall’equazione (2.2), si costruisce il residuo 


eL =r— AZ (2.3) 
L’applicazione delle tecniche ai minimi quadrati — Least Square, LS — o ai minimi quadrati pesati 


— Weighted Least Square, WLS — all’equazione (2.3) porta ai metodi Linear Least Square (LLS) e 
Weighted Linear Least Square (WLLS). 


2.1 Approcci Lineari 


L’idea alla base degli approcci lineari è quella di convertire il modello di osservazione delle varie 
strategie viste nel precedente capitolo in un set di equazioni lineari affette da un rumore a media nulla e 
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di entitá modesta rispetto alla grandezza misurata, cioé caratterizzato da una varianza sufficientemente 
piccola. Una volta che il metodo di localizzazione in esame viene convertito in una forma equivalente 
lineare, la corrispondente funzione di costo risulta sicuramente unimodale per definizione di residuo 
di un sistema lineare (sia esso sovradimensionato o meno), pertanto una soluzione globale del sistema 
€ sempre garantita, nonché conseguentemente una localizzazione verosimile. 

Nel prosieguo, verranno trattato due approcci lineari, ovvero l’approccio lineare ai minimi quadrati, 
o Linear Least Square (LLS), [26, 27] e approccio pesato corrispondente, ovvero il Weighted Linear 
Least Square (WLLS), [28-31] 


2.1.1 Linear Least Square - LLS 


Nell’ambito degli approcci Linear Least Square, il modello di osservazione viene trasformato in un 
sistema lineare sovradimensionato e quindi risolto con le ordinarie tecniche dei minimi quadrati. Nel 
fare ciò, ci si serve di variabili intermedie funzioni della posizione incognita. I passaggi necessari per 
le specifiche implementazioni dei vari approcci sono sviluppati per i singoli metodi di localizzazione, 
pertanto non è possibile generalizzare la strategia; sarà invece necessario operare singolarmente su 
ogni specifico approccio. 


Trilaterazione 


Per convertire la trilaterazione delle distanze in una approccio LLS, si opera generalizzando quanto 
descritto nel caso ideale privo di rumore nell’equazione (1.3) e seguenti, andando a linearizzare in 
modo diretto l'eguaglianza quadratica basandosi sul teorema della corda. Si riparte dell’equazione 
(1.9), valida per il generico sensore l e qui riportata per comodità: 


TRI =d+nr = Va x1)? | (y vi)? FNR, En PEN), (2.4) 


si procede quindi elevando a quadrato entrambi i membri, ottenendo: 


Tha = (8-2) + (yy) +2nriy (e 21)? + (yy)? tng, 1=1,2,...,L. (2.5) 


Si definisce a questo punto una nuova grandezza ausiliaria, il contributo del rumore modificato dal 
processo di linearizzazione: 


MRI =n} + 2nrıy (a — 21)? + (y — yi)?. (2.6) 


Utilizzando quest'ultima equazione puo essere scritto un nuovo modello di osservazione che lega 
posizione e stima della distanza rispetto al sensore l: 


rh, = (2 - ay)? +(y- yi)” FMRI- (2.7) 


N 


E conveniente sviluppare ulteriormente il ragionamento, introducendo la variabile di localizzazione 
ausiliare R definita nel seguente modo 
R= +y. (2.8) 


Di conseguenza l'equazione (2.7) diventa: 
rhi =g? + r? — 2212 + y + Yi — 2y/T + MR (2.9) 


che e utile riorganizzare nel seguente: 


2218 — 2yuy+R+mri=th1-% Y, 1=1,2,...,L (2.10) 
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Questo nuovo modello di osservazione, valido per il sensore /-simo, rappresenta già la generica equa- 
zione di un sistema lineare equivalente a quello di partenza, ma stavolta definito in tre variabili: le due 
coordinate posizionali x, y, ed il loro valore assoluto rispetto al centro del riferimento, ovvero R. Ciò si 
rende più evidente andando a definire quattro grandezze in forma matriciale/vettoriale, in particolare 
la matrice A, legata alla costellazione delle ancore: 


—2zı, -2yı, 1 
— 279, —2y2, 1 


A= : (2.11) 
= Wr, 1 
il vettore del rumore linearizzato, composto da L espressioni introdotte dall’equazione (2.6): 
q= mama... Mp]; (2.12) 
ed infine, la matrice B dei termini forzanti, costituita di due parti: 
b = bı — ba, (2.13) 


dove la prima matrice bı viene determinata in tempo reale, durante la fase di localizzazione vera e 
propria, e rappresenta l'effettivo termine forzante: 


rhi 
TR 

E S (2.14) 
TRL 


mentre la seconda matrice bə può essere costruita in una fase iniziale offline e contestualmente alla 
matrice A, prima di eseguire le misure necessarie a determinare la matrice bi: 


La T Ya 
by = (2.15) 
z} +47 


ed infine il gia citato vettore ausiliario delle incognite 


Si è quindi pervenuti a un sistema lineare affetto da rumore che può essere formalmente schematizzato 
nella seguente scrittura: 
Ap+q=b. (2.17) 


Operativamente, per risolvere il sistema si osserva che, nel caso il rumore di misura sia sufficientemente 
piccolo, si può procedere con la seguente approssimazione: 


2nr 1y (x — 11)? + (y - y)? 
2 2 y (x — 22)? + (y — ya)? 


gr ; (2.18) 


2nrLV (a — ESSE +y- ve)? 
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il vettore q, puó pertanto essere considerato a media nulla, ovvero E [q] = 0, per le ipotesi di partenza 
su nr. Alla luce di quest’ultima considerazione, è possibile approssimare l’equazione (2.17) in valore 
medio, nel seguente modo: 

Ap ~ b. (2.19) 


In questa forma, il sistema della trilaterazione puó essere risolto con le normali metodologie dei 
minimi quadrati. Operativamente si definisce il funzionale linearizzato ai minimi quadrati, indicato 
con JLLS,R(P), nel seguente modo: 


Jırs,r(P) = (AB-b) (AP — b) = P'A' Ap—2p' A'b+b'b; (2.20) 


tale funzionale & una forma quadratica in p, il che assicura l’esistenza di un minimo assoluto. Formal- 
mente, il minimo di questa forma quadratica puó essere determinato eseguendo la derivata rispetto 
alla variabile p calcolata in p: 


sr) q (2.21) 
Op P=P 
da cui discende che 
2A' Ap—2A'b=0=> A'Ap-2A4'b=0 (2.22) 
e quindi, in definitiva 
p=(A A) Ab. (2.23) 


In maniera formale la stima posizionale si ottiene quindi dalla prima e seconda posizione del vettore 
ausiliare p, ovvero: 
2 = [2,9] = [P1,p2]" (2.24) 

Questo modo di procedere è però sub-ottimo, poiché ignorare p3 dopo che esso è comunque stato 
calcolato risulta uno spreco a livello informativo. L'approccio pesato affrontato nella prossima sezione 
fornisce l’opportunità di sfruttare anche questa quantità. 

In alternativa, esiste un secondo modo per linearizzare il modello di osservazione della triango- 
lazione che consiste nell’eliminare R alla radice semplicemente sottraendo due equazioni del modello 
(2.10). Senza perdere di generalità si prenda la prima equazione come riferimento e la si sottragga da 
ogni successiva, ovvero 


22,2 — 2yy + R+ mp1 = TR — xi — y (2.25a) 


2112 — 2yy+R+mr1= rha - gi —y (2.25b) 


ottenendo quindi (L — 1) equazioni linearmente indipendenti del tipo: 


2(xı — x1)£ — 2(yı — yı)y+ mri — MRI = 
2 2 2 2 2 2 = 
TREN TRI tTI ty, 1=2,3,...,L. (2.26) 


Queste possono essere riformulate in maniera matriciale nel seguente: 
Ax+q=b (2.27) 


in questo contesto, le matrici A, q e b sono riformulate nella seguente maniera: 


A (2.28) 
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il termine di rumore 
q = [MR 2 — MR, MR3 — MRI; --.MRLTMRA; (2.29) 


e il termine forzante, ancora una volta scomponibile in contributi variabili con la misura e contributi 
fissi dipendenti dalla sola costellazione: 


2 2 2 2 2 2 2 2 2 2 2 2 
TR2 7 %2— Y2 TRırtTtıry TR2 7 TRA (x3 + ya) — (xi +: 97) 
2 2 2 3 2 2 2 2 2 2 2 2 
TR,3 — %3 — Y3 — TR, T T1 T Y1 TR,3 TRA (123 + y3) - (21 +yi) 
b= . = . = y . (2.30) 
2 2 2 2 2 2 2 2 2 2 2 2 
PRLTYL YLT TRATA Y TRLE-TRA (ei + yp) — (21 +yí) 


Pertanto, assumendo che il rumore q sia sufficientemente piccolo, si procede come nel caso precedente, 
andando ad ottenere cosi la stima posizionale sfruttando le proprieta delle forme quadratiche: 


T=(A'A)!A'b (2.31) 


che stavolta risulta direttamente coincidente con la stima cercata, a differenza di quanto avvenuto nel 
caso precedente. 

Quanto detto in questa sezione vale in maniera equivalente per la trilaterazione basata sui tempi 
di volo, che stima il range d; a partire dai ToF, come descritto nel precedente capitolo. Pertanto a 
livello algoritmico quanto affermato fin qui vale in tutto e per tutto, applicando la semplice operazione 
di moltiplicazione per la velocità della luce laddove necessario. 


Trilaterazione differenziale 


In questo caso si riparte dall’equazione (1.32), riscritta isolando nei due membri i fattori contenenti 
la radice quadrata: 


Tora + y(x- 2%)? + (y - y1)? = Y (2-20? + (y— y)? +npRı (2.32) 


si definisce il vettore di rumore modificato: 


mpri = br + 2npr IV (a — 21)? + (y — yo), (2.33) 


e si introduce inoltre la variabile ausiliaria R, che stavolta assume la forma 


Ri = y (a — 31)? + (y — y1)’. (2.34) 


Elevando a quadrato entrambi i membri si ottiene 


2(x1 — xı) (£ — 21) + 2(y1 — Y) y — y1) + 2rpR Ri + MDR = 
rhbri (£1 — 11) -M-y) 1=2,3,...,L. (2.35) 


Si definisce quindi la matrice A nel seguente modo: 


Tt1=Y YY —TDR2 


Tı — t3 Yo Y3 —TDR3 (2 36) 


11 = IL Y1—YL —TDRL 


quindi il vettore dei termini noti 


p= [£ — z1, y — Yi, Ri]!, (2.37) 
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il vettore del rumore 
q=[mpr2,Mprs)---,MpDRI]', (2.38) 


e il termine forzante: 


rbro— (£1 — 22)” — (y1 — ya)? 


rT (1 — 23)? — (y — ya)” 
ba | i 1 (2.39) 


ron (een) ==) 


Come in precedenza, nell’ipotesi in cui il termine di rumore q sia sufficientemente piccolo, la soluzione 
del sistema risulta: 
p=(A'A)!A'b (2.40) 


da cui, in definitiva: 


E = [2,9] = [Pi + z1, P2 + y1] - (2.41) 


Come nota finale, vale la stessa considerazione fatta al termine della sezione precedente: la stima 
basata sulla differenza dei tempo di volo, ovvero l'approccio noto il letteratura come TDoA, viene 
implementato esattamente come descritto in questa sezione, a patto di scalare il tutto con il termine 
c della velocità della luce. 


Triangolazione 


L’approccio localizzatorio basato sulla triangolazione di rette differisce in maniera abbastanza netta 
da quello basato sulla triangolazione delle distanze e per molti versi ne risulta il duale. In questo caso, 
si riparte dalla relazione generica che lega il valore degli angoli visti dalle ancore alla posizione del 
target rispetto alle coordinate delle stesse: 


raı = Q tnaı = tan”? (=) + NA; (2.42) 
TL XI 


riordinando i due membri si ottiene facilmente la seguente uguaglianza: 


sin (ras — nat) = Y — Y 


ta — = 2.43 
n (rad = nan) cos (rai nal) UL] Ge) 

che puo essere riorganizzata nel seguente modo: 
sin (raı -naı)l® — 21) = cos (rar -naı)ly — ya). (2.44) 


Nell’ipotesi di rumore sufficientemente contenuto, si può senza dubbio ipotizzare che sin (na ı) = 
na, e cos(na:) = 1. Sfruttando queste due approssimazioni, si può semplificare notevolmente il 
modello espresso dall'equazione (2.44) attraverso le formule di somma e differenza delle funzioni 
trigonometriche: 


sin (Tar — na) = sin (ra) cos (4,1) — cos (ra,ı) sin (na) = sin (ra 1) — cos (ra na. (2.45a) 


cos (ra, — nai) = cos (ra) cos (na) + cos (74,1) cos (na) = cos (ra) + sin (ra ı)naı- (2.45b) 


Sostituendo le equazioni (2.45a) e (2.45b) nel modello espresso dall’equazione (2.44) si perviene al 
modello linearizzato: 


sin (ra Ju — cos (ra 1)y + ma, = sin (ra,ı)zı — cos (ra) yi, (2.46) 
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dove é stato utilizzata la seguente definizione del contributo di rumore portato in forma lineare: 


mal = na. |(&ı — x) cos (ra 1) + (yı — y) sin (ra 1) (2.47) 


L’equazione (2.46) può essere riformulata come un sistema lineare attraverso le definizione della 
matrice A: 


sin (rai), —cos(ra1) 
sin(r42), —cos(ra.z) 

A= l i (2.48) 
sin (rar), —cos(ra.c) 


del vettore termine forzante b: 


sin (ra,1)11 — cos (rA,ı) Y1 


sin (74 2)12 — cos (ra 2)ya 
pe (2.49) 


sin (ra.L)ci — COS (TA L)YL 
e del rumore linearizzato q: 
F 
q= [ma MA, 23+- ,Ma,L] . (2.50) 


Pertanto, come nel caso della trilaterazione, si può scrivere il sistema linearizzato nella forma 
Ax+q=b (2.51) 


per il quale si può affermare che il rumore è a media nulla dal momento che E [q] = 0. Sulla base 
di questa considerazione si può procedere come nel caso precedente e identificare la soluzione con la 
seguente inversione formale del sistema lineare: 


T=(A'A)!A'b (2.52) 


2.1.2 Weighted Linear Least Square - WLLS 


Gli approcci di tipo LLS sono molto semplici da implementare, ma pagano tale apprezzabile qualità 
con prestazioni sub-ottime. Infatti, soltanto nel caso in cui i disturbi di misura siano indipendenti e 
identicamente distribuiti i metodi LLS raggiungono l’ottimo prestazionale. Nella realtà invece l’effetto 
del rumore è diverso sulle varie stime, il che significa che, a priori, è possibile attribuire un certo livello 
di affidabilità alla singola misurazione r;, e quindi sfruttare questa metrica nella ricombinazione delle 
misure. 

L’accuratezza dei metodi LLS può essere migliorata includendo una matrice di peso, indicata con 
W. A partire dalla generica equazione linearizzata affetta da errore in forma lineare, ovvero 


Ap+q=b (2.53) 
si definisce il funzionale quadratico pesato con la seguente forma, indicata come Jw yz Ls(0): 
JwLLs(p) = (Ap— b)W(Ap — b). (2.54) 


In accordo con la teoria del miglior estimatore lineare non polarizzato (Best Linear Unbiased Estimator, 
BLUE), [28], la matrice W deve essere scelta uguale all'inverso della matrice di covarianza di q, in 
modo da ottenere la prestazione asintoticamente ottima. Qualitativamente, la matrice W introduce 
l’effetto della statistica del rumore nella stima di r e permette di migliorare l’algoritmo di localizzazione 
aggiungendo delle informazioni sull’affidabilità della singola equazione del sistema lineare. Data la 
natura della matrice W, non è possibile generalizzare l'approccio pesato, che si specializza a partire 
dalla specifica strategia di linearizzazione applicata allo modello di osservazione selezionato. 
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Trilaterazione 


Nel caso della trilaterazione, considerando l’espressione del rumore per l’equazione linearizzata 
> P P q > 


OVVero 
2nrıy (a — 21)? + (y — m1)? 
Anray (x — 22)? + (y — ya)? 


q= (2.55) 
Mp, y (a — 21)? + (y — yL)? 
si può costruire la matrice W nel seguente modo: 
1 a -1 
W = [E[ga']] * = [diag (40% 47,402. ad, ... 40h rda) - (2.56) 


Dal momento che il vettore dei valori d; non è ovviamente noto in forma esatta, la migliore strategia 
per la costruzione della matrice W passa attraverso la sostituzione di questi ultimi con la loro migliore 
stima disponibile, ovvero r;. Pertanto: 


? 2 a 2 2 
oRırı 407 213 40h ri 


1 1 1 
w = aie (73 5 ies ) (2.57) 


Si può da questo punto procedere con l’approccio classico, determinando esplicitamente il minimo del 
funzionale quadratico nella sua variante pesata come previsto dall’equazione (2.54), ovvero: 


p= argmin (JwLLs,r(p)) = (A'WA) 4 Wb. (2.58) 
P 


Come nel corrispondente caso LLS, la stima posizionale si ottiene quindi dalla prima e seconda 
posizione del vettore ausiliare p, e cioè: 


T= Ec, y) = Pi. Dal". (2.59) 


Questo approccio fornisce risultati migliori di quelli determinati dall’equazione (2.23), al prezzo di 
una semplice operazione aggiuntiva di scarso peso computazionale. 


Trilaterazione a due passi 


Il risultato ottenuto utilizzando il pesaggio può essere ulteriormente migliorato sfruttando la stima 
di R = yx? + y?, fornita da p3, disponibile anch’essa a costo nullo. Aggiungendo un secondo stadio 
di elaborazione all’algoritmo, che diviene così un metodo a due passi, si ottiene un risultato ancora 
più raffinato semplicemente applicando alcune operazioni aggiuntive, [29]. 

Si noti infatti che se £ è una stima sufficientemente accurata di x, la differenza di quadrati può 
senza dubbio essere sviluppata nel seguente modo: 


(61) — a? ~ 2a (fi — 2) (2.60a) 


(P2)? — y? ~ 2y (Pa — y). (2.60b) 


Si può quindi procedere col costruire il sistema ausiliario nella variabile z = [x?, y?] ', da determinarsi 
sulla base delle stime ottenute dalla soluzione precedente, che diviene cosi il passo 1 di un sistema in 
due step. Il vettore z è il risultato del sistema lineare espresso da: 


h=Gz+w (2.61) 
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che coinvolge un termine forzante costruito a partire dalla soluzione completa del funzionale (2.58): 
A (a2 a 
h= Mm , (P2) Ps, (2.62) 
e la seguente grandezza ausiliaria: 


1 0 
G=|o 1|, (2.63) 
1 1 


nonché il termine “rumoroso” derivante dalle espressioni (2.60): 
P 


w-22(f-2), 2-4), PR]. (2.64) 


Per determinare la matrice di covarianza del vettore w si utilizza la teoria dell’estimatore BLUE, da 
cui si puö concludere che [32]: 


lo -y Rd 29-37] =(ATWA) ; (2.65) 


pertanto, utilizzando i risultati ottenuti dalle equazioni (2.64) e (2.65), la matrice di peso ottimale 
risulta in questo caso pari a 


$ = [diag(2a, 2y, 1)(ATW A)! diag(2x, 2y, 1)| a (2.66) 
che in definitiva porta all'ottima soluzione pesata 
2 = (Gece) "Glen. (2.67) 
Finalmente, la nuova stima per la posizione risulta 
@ = [sign (Ai) VA. sign) VA] (2.68) 


dove sign(-), ovvero la funzione segno, si rende necessaria dal momento che pj e Py sono i quadrati di 
rey. 


Triangolazione 


In questo caso, la determinazione della matrice di peso sfrutta alcuni risultati validi per la trilate- 
razione. Innanzitutto è necessario ipotizzare di avere a che fare con rumori sufficientemente contenuti 
per poter supporre che 


x — x, = dı cos (di) (2.69a) 
y — yı = di sin (61). (2.69b) 


In queste condizioni il vettore di rumore modificato si può semplificare in 
ma = na, |(zı— x) cos (ra) + (yı — y) sin (r41)] © na |(xı — x) cos (61) + (yı —y)sin(d)] (2.70) 
da cui, sfruttando le due equazioni (2.69): 


MA, = -na,.l-dı COS (di)? = di sin (1)"] = NAL dı. (2.71) 
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La matrice di peso può quindi essere determinata come nel caso della trilaterazione, ovvero: 


| E u tees ) (2.72) 


2 2222 Dia 1222 2 
Ad 042% Andi 


Sfortunatamente, le distanze d; non sono disponibili in questo contesto e non sono neppure legate in 
maniera diretta alle grandezze misurate su cui si basa l’algoritmo come nel caso della trilaterazione, 
dato che in questo approccio sono invece gli angoli di arrivo del segnale a svolgere il ruolo di quantitä 
in ingresso. Pertanto, per costruire le matrici di pesaggio è necessario ricorrere ad una strategia 
alternativa, di tipo iterativo. Questa può essere sviluppata nei seguienti passi: 


1. si determina & utilizzando l’approccio LLS della triangolazione, descritto nel paragrafo 2.1.1; 
2. si stimano le distanze dı utilizzando la la stima £ del punto precedente; 

3. si costruisce la matrice W utilizzando le stime d; e le statistiche delle stime angolari 04,1; 

4 


. si determina una nuova stima di £ utilizzando l’approccio pesato; 


or 


si ripetono i passi 2-4 fino a che la stima della posizione non è sufficientemente raffinata. 


Anche in questo caso si ottengono risultati sensibilmente migliori rispetto che con la versione non 
pesata, a prezzo di una piccola complicazione computazionale. 


2.2 Esempi di triangolazione basata su approcci lineari 


In questa sezione sono riportati alcuni test di localizzazione basati sugli algoritmi presentati. Si 
fa riferimento in generale ad uno dominio di 10 metri per 10 metri, monitorato da 4 ancore poste ai 
quattro vertici del quadrato, come rappresentato in Figura 2.1. In questo dominio il target è libero di 
muoversi in qualunque punto. 

Con riferimento alla Figura 2.1, la posizione del target T viene stimata sulla base delle distanze 
individuate rispetto al set di ancore disponibili. Per fissare le idee il target è inizialmente posizionato 
nel punto x = [3,6], e la stima delle distanze è stata modellata come affetta da un rumore bianco 
gaussiano di varianza 07. Questa varianza è supposta proporzionale al valore del range dı, sulla base 
di considerazioni che trovano amplio riscontro nella pratica. Per valutare un ampio spettro di casi, si 
è reso variabile il rapporto segnale rumore (SNR), definito in modo da soddisfare la proporzionalità 
con la distanza, ovvero: 

di 
> 


0; 


SNR= (2.73) 
Avendo a disposizione questa libertà si è simulato un valore di SNR non costante nel dominio, ma 
che invece migliora in vicinanza delle ancore, come di fatto succede in un contesto reale. Per conferire 
a questo test un valore statistico, la stima posizionale è stata eseguita 1000 volte, e in ognuna di 
queste realizzazioni simulate il rumore bianco è stato generato casualmente ma sempre soddisfacendo 
la relazione di proporzionalità e mantenendo costante il rapporto segnale rumore in ogni prova. Infine, 
vari test sono stati condotti variando il valore di SNR a partire da valori estremamente bassi, per 
modellizzare situazioni critiche, fino a valori che rasentano l’idealità. 

Per ognuno dei casi analizzati la posizione del target è stata stimata utilizzando l’algoritmo LLS 
di tipo 1, l’algoritmo WLLS semplice e la variante WLLS a due passi. In Figura 2.2 si può apprezzare 
come, a parità di rapporto segnale-rumore, l’algoritmo caratterizzato dalle le migliori prestazioni sia 
sempre l’algoritmo WLLS a due passi, mentre il peggiore l’algoritmo LLS base. È da rimarcare come 
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Figura 2.1: Dominio di localizzazione di 100 m? e relativa costellazione di 4 ancore, posizionate ai vertici 
di un quadrato di lato 10m. Due target sono stati presi un considerazione: il primo in posizione [3,6], e il 
secondo in [2,1]. 


la variante WLLS a singolo passo sia praticamente sovrapposto a quest’ultimo per valori di SNR 
maggiori di 10 dB. 

Un secondo test è stato condotto considerando la posizione [2, 1], e i relativi risultati sono presentati 
in Figura 2.3. In questa seconda prova i risultati sono peggiori in termini di accuratezza rispetto a 
quelli ottenuti in precedenza con lo stesso valore di SNR. Inoltre, appare più evidente la differenza tra 
i vari algoritmi, con LLS che si conferma il peggiore mentre WLLS a due passi che fornisce risultati 
sensibilmente migliori di WLLS singolo. Asintoticamente, tra il migliore e il peggiore risultato ci sono 
circa 5 dB di differenza. Questa situazione non è casuale, ma corrisponde al fatto che la posizione 2 è 
più critica della posizione 1, per motivi che saranno resi più chiari nel prossimo capitolo. 


2.3 Approcci non lineari 


Nell'ambito degli approcci non lineari si cerca di determinare la posizione del target minimizzando 
direttamente le vari equazioni affrontate nel capitolo precedente, viste come funzionali da ottimizzare. 
La convergenza globale di questi approcci non è pero garantita a priori, in quanto le relative funzioni 
di costo sono in generale multimodali, ovvero possono presentare dei valori corrispondenti alla mini- 
mizzazione del residuo che però non corrispondono al valore effettivo della posizione del target; ciò 
accade a causa della presenza del rumore di misura. 

Paradossalmente, i metodi nonlineari sono comunque più semplici da implementare rispetto a quelli 
lineari, perché non richiedono la conoscenza specifica dei meccanismi con cui la stima della posizione 
viene determinata, cosa invece necessaria per ottenere una soddisfacente linearizzazione dell’equazione 
del modello di osservazione. 

Il metodo ai minimi quadrati di tipo non lineare, ovvero Non-Linear Least Square (NLS), in 
particolare, è il più semplice da implementare ed è una scelta pratica nei casi in cui non sia disponibile 
alcuna informazione sulla statistica del rumore. L’implementazione del metodo NLS infatti sfrutta 
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Figura 2.2: Stima della coordinate cartesiane del target 1 in Figura 2.1 mediante lateralizzazione imple- 
mentata con gli algoritmi lineari descritti nei paragrafo 2.1.1 e 2.1.2: LLS, WLLS a singolo passo e WLLS 
a due passi. 


strategie che sono molto generali, andando ad applicare risultati della teoria dell’ottimizzazione nella 
loro forma pit teorica. 

D'altra parte, il metodo basato sullo stimatore a massima verosimiglianza, ovvero Maximum Like- 
lyhood (ML) utilizza in modo costruttivo la matrice di covarianza del rumore, ed è ottimale nel senso 
che la sua prestazione può eguagliare il limite teorico sancito dalla teoria dell’informazione, come sarà 
chiarito in seguito. Può inoltre essere considerato la versione pesata del metodo NLS, nello stesso 
modo in cui i metodi lineari di tipo WLLS estendono quelli LLS utilizzando la teoria del migliore 
stimatore linearizzato. 


2.3.1 Nonlinear Least Square - NLS 


Il metodo NLS minimizza la funzione di costo di tipo ai minimi quadrati partendo dalle equazioni 
(1.7), (1.41) e le simili implementazioni dei modelli di osservazione. 

A titolo di esempio si analizza il caso della trilaterazione dei range dı, ma l'estensione alla trian- 
golazione e tutti gli altri approcci possibili è estremamente diretta. Basandosi sull’equazione (1.7) la 
funzione di costo NLS, indicata con Jnzs,r(&£), può essere scritta nel seguente modo: 


InLs,r(®) = (rr - fr(@))' (rr - fr(@)) (2.74) 
o alternativamente, esplicitando la sommatoria delle componenti dei vettori: 
2 


InLs,R(È) = 5 (rra =e =a F= u”) ; (2.75) 


l=1 


per cui il problema della localizzazione si traduce formalmente nel determinare la posizione 2 cui 
corrisponde il più piccolo valore del funzionale Jnrs,r. Dal momento che il funzionale è sempre 
definito positivo, è formalmente lecito scrivere che: 


T = argmin (Jnzs,R(&)) (2.76) 
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Figura 2.3: Stima della coordinate cartesiane del target 2 in Figura 2.1 mediante lateralizzazione imple- 
mentata con gli algoritmi lineari descritti nei paragrafo 2.1.1 e 2.1.2: LLS, WLLS a singolo passo e WLLS 
a due passi. 


Determinare è non è un compito semplice in quanto nel dominio JyLs,r vi sono minimi locali oltre 
al minimo assoluto. Fondamentalmente, ci sono due strategie per risolvere (2.76). La prima tenta di 
eseguire un’esplorazione globale utilizzando la ricerca della griglia o ricerche casuali, come algoritmi 
genetici [33] e ottimizzazione a sciame di particelle, noto con il nome di particle swarm nella letteratura 
internazionale [34]. 


L’altra strategia possibile corrisponde ad una ricerca locale, ovvero un algoritmo iterativo basato 
su una stima iniziale della posizione, indicata con 2, dove 0 si riferisce alla stima iniziale, anche 
grezza. Dal momento che tutti gli algoritmi qui trattati si riferiscono al caso indoor, una stima iniziale 
puo corrispondere al centro del dominio di localizzazione, o eventualmente alla posizione di una delle 
ancore scelta casualmente, o infine al risultato il di un altro metodo di localizzazione, ad esempio uno 
dei metodi lineari descritti nella sezione precedente. In quest’ultimo caso è consigliato implementare 
i metodi LLS, computazionalmente semplici, piuttosto che i più sofisticati ma più complessi metodi 
WLLS. Se questa stima iniziale 2" risulta sufficientemente vicino a x, si suppone che la soluzione 
esatta x possa essere determinata in maniera abbastanza rapida con una procedura iterativa. 

Qui affronteremo i metodi di ricerca locale, in particolare i metodi Gauss-Newton, Newton-Raphson, 


e il metodo di massima discesa, noto internazionalmente con il nome di Steepest Descent. Per tecniche 
di ricerca locale più avanzate, si rimanda a [35]. 


Gauss-Newton 


Il metodo Gauss-Newton sfrutta una considerazione molto semplice per costruire la successio- 
ne delle stime. Per risolvere il problema espresso dalla (2.76) si suppone in questo contesto che la 
stima iniziale ° sia così vicina alla soluzione da permettere uno sviluppo di Taylor nell’intorno di 
quest’ultima. In questa ipotesi si può pertanto scrivere 


fa(2) = fr(@) + G*@ — 2") (2.77) 
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dove G* è pari a G fr(x)), ovvero la matrice Jacobiana di fr(&°), calcolata in af. Nel caso della 
trilaterazione, questa matrice assume la seguente forma esplicita: 


y le-21)?+W-yı)? y (ea) ?+ (yy 1)? 

dx ’ Oy 

Af (am) +Hy=y2)? Oy le-22)?+(y-y2)? 
G(fa(z)) = di iù 5 


OV(e-c2)?+(y-y2)? 0V(e-c2)?+(y-y2)? 
dx ’ Oy 


T—L1 Yy—T1 
((w@—a1)?+(y—y1)?) 7?’ (ee) He) 
TL XI 


1 
((w@—w2)?-+(y—ya)?) 1/7? (0-82)? +y-Y2)2)"’? 


(2.78) 


TEEL Yy-TL 
(2021 )?.+(y-yL)23)1/2? ((e-21)?+Y-Yr)?)/? 


Cambiando il metodo di localizzazione cambia la funzione di osservazione, e conseguentemente la 
matrice Jacobiana, ma il procedimento resta identico. 

Sostituendo lo sviluppo prodotto nell’equazione (2.77) nel funzionale espresso dall’equazione (2.74) 
si ottiene: 


A 30 NS „0 ki 0 
Inus,n(@) = (ra-(fa@)+6"@-2)) (re-(fr@)+6*@-3))). (279) 
Si definisce quindi un nuovo termine forzante, basato su quello originale e sulle grandezze ottenute 
grazie alla stima iniziale x°: 
rr'=rr-fr(®)+G"3°; (2.80) 
considerando questo termine, l’equazione (2.79) può essere riscritta nel seguente 


rale (re = cè) (re = c'è) l (2.81) 


Questa forma è lineare in £ e corrisponde quindi ad un problema ai minimi quadrati sovradimensionato. 
Come discusso nella sezione precedente, la soluzione canonica del sistema risulta: 


= (GA G*)-1G*' rp’; (2.82) 
tuttavia, utilizzando l’equazione (2.80), questa espressione risulta pari a 
@ = (GFT GEIGE (rr ER) + G'a.) l (2.83) 
e quindi riordinando i termini e sfruttando le proprietà delle matrici: 
= 2° +(G*'G*)-1G*" (rr = fr(@°)) (2.84) 


Osservando questa espressione si evince come a partire dalla stima iniziale & si sia pervenuti ad 
una nuova, basata sull’attuale quantita residuale (rr — fr(&°) , che esplicitamente quantifica la 
distanza della stima disponibile nei confronti dalla soluzione reale. Inoltre l’espressione matriciale 


T I ; RE : : : HR: ; ER 
(GF Er PE fornisce la direzione da seguire per aggiornare questo residuo nella direzione in cui si 
trova la soluzione reale. Si puó dunque pensare di generalizzare il metodo, reiterando il processo sulla 
stima di £, generalizzando la regola di aggiornamento: 


pa (Cr) EUREN)" (re fr) (2.85) 
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Newton-Raphson 


Il metodo di Newton-Raphson generalizza ancor di piú l’approccio Gauss-Newton. In questo caso, 
a partire dalla stima iniziale x, la regola di aggiornamento per la posizione risulta la seguente: 


et) = è HU (Inrsr(®"))V(InLsr(®")), (2.86) 


dove VINLS.R(& ))è il gradiente di Jvrs,r(&£)), valutato al passo k*?: 


OJNLS,R(®) 
V(JInzs,r(®)) = | ar, re) (2.87) 
Oy 
e H(Jnıs.r(&")) è la matrice Hessiana, anch'essa valutata al passo k*”: 
2 
Da) le) Inrs nle), ə? Inter (e) 
H(JNLs,r(2)) = = bada (re E miste (2.88) 
dxdy ’ Oy? 


Nel caso della trilaterazione queste grandezze si specializzano in: 


O InLs.R(®) _ 5 2 (z — a)? (rr. -= y(x = 11)? + (y — y1)? ) (y — y) ; 
( y 


= EACT APRU ((e-2)?+(W- 2) : 
(2.89a) 
Chests el, = y)? (rei - V@— a? FU- yi)? ) (a — 2)? 
ep A RSE 2 (œ= 21)? + (w= 40) Di 
(2.89) 
e inoltre í 
O Inrs,r(®) _ O Ino R(E) _ > Fr i(2— any — ya) os 
N Ordy ala)? y ye 
mentre il vettore gradiente risulta 
A nie) | L (rriz a+) (e w] 
vu dire lg eo ol 
a) fera L (rra (z—21)?+(y vi) DI (2:30) 
El (=a P Hun 


Steepest-Descent 


Infine, viene presentato il procedimento iterativo per il metodo della massima discesa, noto in 
letteratura come metodo Steepest-Descent. Per questo, la formula di aggiornamento della stima risulta: 


et = 2 — uV (Jnrs,R(@*)), (2.91) 


dove u è una costante positiva che controlla la convergenza e le stabilità. In generale, più alto il valore 
di u e più rapida è la convergenza, ma anche più elevato il rischio dell’instaurarsi di oscillazioni; nella 
pratica è necessario scegliere un valore adeguatamente basso per valutare la stabilità dell’algoritmo. 
Quale che sia il metodo iterativo scelto, e quindi la corrispondente equazione di aggiornamento 
(2.85), (2.88) o (2.91) , è comunque necessario stabilire un criterio di arresto. Tipicamente, per questi 


¿eL _ ak N 
metodi si stabilisce un massimo per il raffinamento iterativo, nella forma di ||2"* llo < e, dove e è 
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un numero positivo adeguatamente piccolo. Inoltre, è tipicamente necessario stabilire anche il massimo 
numero di passi, onde evitare lo stallo nel caso di mancata convergenza, allorquando la soluzione oscilli 
tra due stati equidistanti dal valore reale. 

Per trarre un bilancio, il metodo di Newton-Raphson e il metodo di Gauss-Newton sono molto 
rapidi, ma pagano questa rapidità con una grossa onerosità computazionale, dal momento che è im- 
plicata l’inversione di una matrice. Di converso, il metodo Steepest Descent è stabile, ma il suo tasso 
di convergenza è basso. Quest'ultimo metodo può essere considerato una forma di approssimazione 
dell’equazione 2.86, in cui il calcolo della matrice Hessiana è omesso. 


Applicazione ad altri approcci 


Nel caso si applichi un metodo della stima della distanza che sfrutti la potenza ricevuta, quali quelli 
basati sull’inversione della stessa — i.e. Received Signal Strength Indicator, RSSI — è possibile riscrivere 
il tutto sulla base della differente definizione. Formalmente quindi, il nuovo funzionale, indicato come 
INLS,RSSI (2) , risulta 


Inıs.nssı(®) = YO (assi +aln (VE a)? + Gw) = 


l=1" 


(rrsst — frssi(®)! (rrssı - frssı(®)), (2.92) 


pertanto la corrispondente stima NLS risulta formalmente individuata dal seguente: 
T= argmin (InLs.rssi(€)) . (2.93) 
zT 


Lo stesso approccio puo essere facilmente generalizzato a tutti gli altri metodi, anche alla trilaterazione 
differenziale o alla triangolazione. La differenza sta nel calcolo delle rispettive metriche e matrici, ma 
in ogni caso questa operazione puó essere compiuta in forma chiusa, data la struttura delle funzioni 
di osservazione coinvolte. 


2.3.2 Maximum Mikelihood - ML 


Risultati migliori per quanto riguarda la stima posizionale si possono ottenere partendo dalla 
conoscenza della distribuzione dell'errore. Gli approcci basati sulla massima verosimiglianza, ovvero 
quelli raggruppati sotto la dicitura Maximum Likelihood, utilizzano la conoscenza della statistica 
dell’errore per ottenere risultati ottimi, nel senso che l’errore che si commette sulla stima della posizione 
e compatibile con il limite asintotico per questo tipo di problema — si veda il prossimo capitolo a 
riguardo. Nel caso in cui l’errore sia guassiano a media nulla, la massimizzazione delle equazioni 
corrisponde alle versione pesata dell'approccio NLS. 

Operativamente, si prende in considerazione la densitá di probabilitá associata al modello di os- 
servazione, ovvero la forma definita da (1.13) per quanto riguarda la trilaterazione, ripetuta qui per 
comoditä: i 


(2r)? |Cr 


1 2 
p(rr)= na exp ( 5 (TR d)' CH (rr - a)) (2.94) 
Operando in questo contesto, la posizione del target & identificata come quella che massimizza la 
densità di probabilità. Si preferisce però considerare la versione logaritmica di questa densità di 
probabilità, definita generalmente con il nome score nella letteratura internazionale: 


nora) =n (rare) pad) (2.95) 
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Il primo termine del secondo membro di (2.95) è indipendente dalla variabile x, per cui massimizzare 
l'equazione (2.94) equivale a minimizzare il solo secondo termine (ra — d)' CE (rr — d). Definito 
quindi il funzionale 


Imr,r(®)= (rR - fa(2)) CR (rr - fa(2)) = > = (2.96) 
l=1 Rl 
la stima di tipo ML risulta formalmente: 
T = argmin (JmL,R(©)) . (2.97) 


Confrontando le equazioni (2.96) con (2.74) e (2.75), si puó osservare come, in presenza di un rumore 
bianco gaussiano a media nulla, lo stimatore di tipo ML generalizza il metodo NLS introducendo un 
fattore aggiuntivo al singolo termine della sommatoria, ovvero operando un pesaggio del contributo 
informativo dovuto alla stima l- dee 

Intuitivamente, se la varianza 0? R ı relativa alla misura dell’ancora l-sima è elevata, significa che 
questa misura è da considerarsi poco affidabile, certamente meno affidabile di quella fornita da un 
altro elemento che presenti o più bassa. Da ciò consegue come nella combinazione di tutti i termini 
il contributo delle due ancore debba essere pesato con un valore che tenga conto della loro reciproca 
relazione di forza nella sommatoria, in modo da rendere conto della rispettiva attendibilità. Una misura 
caratterizzata da un rumore molto modesto sarà infatti più affidabile, e quindi il suo peso 1/ oR, 1 sara 
un numero proporzionalmente grande rispetto a quello delle misurazioni inaffidabili. 

Se tutte le stime sono caratterizzate dalla stessa varianza, ovvero se la matrice Cig? risulta pro- 
porzionale alla matrice identità attraverso un mero coefficiente scalare, l'approccio ML si riduce a 
quello NLS generico. Il legame tra ML e NLS è quindi lo stesso che sussiste tra WLLS e LLS, ovvero 
Vapproccio ML è in grado di utilizzare con successo le informazioni aggiuntive che sono implicitamente 
contenute nella statistica del rumore sovrapposto alla misura. 

Operativamente, per determinare la soluzione è necessario ricorrere a un metodo risolutivo non- 
lineare, come quelli affrontati precedentemente. Utilizzando il metodo di Newton-Raphson: 


et) — è" H! (Jur, R(@®))V (JmL,R(@®)), (2.98) 
dove 92 È 
_ JML, R(E 

H(Jut,r(@)) = oT’ (2.99) 


ovvero, nel caso della trilaterazione si ha: 


e (x - 21) (rei - VE- a? + Wu?) (y — u)? 
[H(Jur,n(2) "Za (z — x1)? + (y— y? (oca EN. 

(2.100a) 

3 (y = u)? (rn. y (a = a1)? + (y y) (0 a)’ \ | 

[IH (Jmı,r(&))]2 2=2, oz, \(e-a)?+(y- wm? Ge (y- 2)? ; 

(2.100b) 
inoltre 

= 2rru(e = 21) (y — y) 
[IH (Juı,r(&))ı,2 = [IH (Jmı,r(&))l2,ı = 5 i (2.100c) 


l=1 Thi (=P. + (y — yyy? 
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e infine 


(rra (121)2+(y vu) (e x) 
or, (@-21)?+@y-91)?)°7/? 


F agli 2.101 
(Jmz,r(&)) > (ra. (2-21)2+(y vi?) w vi) 
= oF, ((@—21)? +(y—)?) >? 


E inoltre possibile procedere con gli approcci Gauss-Newton, sfruttando una forma che generalizza 
quella vista in (2.80): 


2 = F" + (G' (fa(2")C 7 G(fa(2") G! (fr(@"))CR (re - fr(@")): (2.102) 
e infine, mediante il procedimento iterativo per il metodo della massima discesa: 
gt! = è — uV(JImr,r(&"), (2.103) 


in cui valgono le stesse considerazioni viste precedentemente sul valore di u e sui criteri di arresto. 


2.4 Esempi di triangolazione basata su approcci non lineari 


Anche per il caso degli algoritmi non lineari è opportuno fare un confronto. Innanzitutto è necessario 
valutare come, fissate le condizioni al contorno (ovvero il rapporto segnale-rumore), i tre algoritmi 
iterativi si comportino rispetto alla rapidità di esecuzione. A questo proposito, con riferimento allo 
stesso setup modellato in Figura 2.1, la figura Figura 2.4 mostra il risultato dei tre algoritmi Gauss- 
Newton, Newton-Raphson e Steepest Descent, quando il target si trova nella posizione x = [3,6], 
ovvero la stessa posizione su cui sono stati esegui i test con gli approcci lineari nella sezione 2.2. 
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Figura 2.4: Stima della coordinate cartesiane del target utilizzando i tre algoritmi NLS descritti nel 
paragrafo 2.3.1: Gauss-Newton, Newton-Raphson e Steepest Descent. 


Per tutti e tre i casi è stato scelto come punto di partenza il centro del dominio, ovvero il punto 
To = [5,5], per mettere in evidenza la rapidità di convergenza dei tre metodi. In figura Figura 2.4 si 
può apprezzare l'andamento della convergenza dei tre algoritmi per la coordinata x, a sinistra e la 
coordinata y, a destra. Per fissare le idee è stato scelto un elevato livello di SNR, pari a 30 dB. 

Figura 2.5 mostra invece l'andamento degli stessi tre metodi non lineari, ma stavolta implementati 
secondo la filosofia della massima verosimiglianza basata sulla statistica dell’errore, come descritto 
nella sezione 2.3.2. Anche in questo caso si conferma la stessa efficacia dei tre approcci, nel senso 
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che ancora una volta i metodi Gauss-Newton e Newton-Raphson forniscono la stima definitiva in 
un numero ristretto di passi, mentre il metodo della massima discesa impiega almeno una decina di 
iterazioni. 


I I 6.5 T T T I I 

—e— ML: Gauss-Netwon —e— ML: Gauss-Netwon 

—— ML: Newton-Raphson —— ML: Newton-Raphson 
45 —¢— ML: steepest descent |] —*— ML: steepest descent 
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Figura 2.5: Stima della coordinate cartesiane del target utilizzando i tre algoritmi NLS descritti nel 
paragrafo 2.3.2: Gauss-Newton, Newton-Raphson e Steepest Descent. 


L’errore che si ottiene utilizzando gli approcci NLS é paragonabile a quello che si ottiene utilizzando 
Vimplementazione ML. Con riferimento alla Figura 2.6, anche in questo caso si è valutato l’effetto del 
rapporto segnale rumore nella stima della posizione nelle 2 varianti NLS e ML. In particolare si fa 
riferimento all’implementazione Newton-Raphson, dal momento che, fissato il tipo di pesaggio, i tre 
algoritmi forniscono lo stesso risultato, mentre cid che cambia é il numero di passi necessari per 
ottenerlo. Il comportamento é simile nei due casi, producendo curve praticamente sovrapposte, anche 
se si può apprezzare come quella riferita al metodo ML sia marginalmente migliore. 

Ripetendo l’intera analisi per una diversa posizione del target, si osservano alcune differenze. 
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Figura 2.6: Confronto tra tre le implementazioni Newton-Raphson di tipo NLS, descritto nel paragrafo 
2.3.1 e di tipo ML, descritto nel paragrafo 2.3.2. 
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Considerando di nuovo il punto « = [2,1], come fatto nel paragrafo 2.2, i risultati degli approcci 
NLS sono riportati in Figura 2.7, mentre quelli per gli approcci ML sono riportati in Figura 2.8. 
Si può osservare come la stima della posizione sia più difficoltosa per tutti gli algoritmi, in tutte le 
sei combinazioni di NLS/ML e GN/NR/SD. In questo caso però emerge che i metodi ML sono in 
vantaggio rispetto a quelli NLS, almeno per bassi valori di SNR, come si può apprezzare osservando 
Figura 2.9, che mostra l’errore quadratico medio in funzione del rapporto segnale rumore. Questa 
differenza dipende dal migliore utilizzo che gli algoritmi ML fanno delle informazioni disponibili. 
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Figura 2.7: Stima della coordinate cartesiane del target utilizzando i tre algoritmi NLS descritti nel 
paragrafo 2.3.1 quando il target si trova in posizione x = [2, 1]. 
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Figura 2.8: Stima della coordinate cartesiane del target utilizzando i tre algoritmi ML descritti nel paragrafo 
2.3.2 quando il target si trova in posizione x = [2,1]. 


2.5 Confronto tra approcci lineari e non lineari 


Per concludere questo capitolo, tutti gli algoritmi presentati sono messi a confronto. In particolare 
Figura 2.10 li mette a confronto quando la posizione del target è x = (3,6). Si può osservare come, 
asintoticamente, tutti i metodi presentino risultati simili e soltanto per bassissimi valori di SNR 
algoritmo WLLS a 2 passi risulti il migliore. Asintoticamente, l’errore quadratico medio è compreso 
entro 1 dB. 
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Figura 2.9: Confronto tra tre le implementazioni Newton-Raphson di tipo NLS, descritto nel paragrafo 
2.3.1 e di tipo ML, descritto nel paragrafo 2.3.2 quando il target si trova nella posizione x = [2, 1]. 


Considerando invece Figura 2.11, i vari algoritmi si distribuiscono in una dinamica di circa 4dB. 
Si può osservare come, asintoticamente, il metodo WLS a 2 passi e il metodo Maximum Likelyhood 
siano praticamente identici producendo il piú basso errore quadratico a paritá di SNR, mentre l'imple- 
mentazione LLS, la pit semplice tra tutte quelle presentate, si conferma come quella con le peggiori 
prestazioni. Le due posizioni æ = [3,6] e æ = [2,1] sono state scelte proprio perché molto estreme 
nella dinamica delle possibilità. Nel prossimo capitolo verrà chiarito il perché di questa differenza di 
risultati. 

Come nota finale a questo capitolo, la medesima successione di prestazioni tra gli algoritmi vale 
anche per le altre implementazioni, in particolare per la triangolazione basata sull’intersezione di rette. 
Cambia eventualmente l’effetto che la posizione del target ha sull’errore di localizzazione osservato, nel 
senso che, fissata la costellazione, la propagazione dell’errore di misura è diversa per i vari algoritmi. 

Per concludere questo capitolo, la tabella 2.1 mostra un confronto tra vantaggi e svantaggi dei 
vari metodi, lineari e non lineari, discussi in questo capitolo. In conclusione, quelli che considerano la 
statistica dell’errore ottengono generalmente i risultati migliori, ma non sempre quest’ultima è nota: 
introdurne una sbagliata può peraltro peggiorare molto il risultato finale. I metodi lineari determinano 
sempre la soluzione globale, ma la loro accuratezza può essere scarsa; d’altro canto, quelli non lineari 
sono generalmente più prestanti, ma anche più onerosi. 


38 INTRODUZIONE AI SISTEMI DI LOCALIZZAZIONE INDOOR 


50 r r r — 1 z 
} —e_ LLS 
4 DR A rn Rt WLLS LÌ 
EN ARTO IO ee —- WLS a 2 passi I 
IN —4— NLS Newton-Raphson 
DOVE as AA E —*— ML Newton-Raphson 


Errore quadratico medio (dB) 
o 


-10 0 10 20 30 40 50 60 
SNR (dB) 


Figura 2.10: Confronto tra tutti i metodi di triangolazione presentati in questo capitolo quando il target 
si trova in posizione æ = [3,6]. LLS introdotto in 2.1.1, WLLS e WLLS a 2 passi introdotti in 2.1.2, 
Newton-Raphson di tipo NLS, introdotto in 2.3.1 e Newton-Raphson di tipo ML, introdotto in 2.3.2. 
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Figura 2.11: Confronto tra tutti i metodi di triangolazione presentati in questo capitolo quando il target 


si trova in posizione @ = [2,1]. LLS introdotto in 2.1.1, WLLS e WLLS a 2 passi introdotti in 2.1.2, 
Newton-Raphson di tipo NLS, introdotto in 2.3.1 e Newton-Raphson di tipo ML, introdotto in 2.3.2. 
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Tabella 2.1: Confronto tra i vari approcci presentati in questo capitolo 


Metodo Vantaggi Svantaggi 

LLS Implementazione semplice Bassa accuratezza 
Soluzione globale garantita 
Nessuna conoscenza pregressa 
richiesta 

WLLS Soluzione globale garantita Richiesta conoscenza delle stati- 
Media/elevata accuratezza stiche dell’errore 

Necessitä di iterazione 

NLS Media/elevata accuratezza Soluzione globale non garantita 
Nessuna conoscenza pregressa Complessità media/alta 
richiesta Necessità di iterazione 

ML Elevata accuratezza Soluzione globale non garantita 


Complessità media/alta 
Necessità di iterazione 
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3. DILUIZIONE DELLA PRECISIONE 


Nel capitolo precedente si sono mostrati vari metodi per la determinazione della posizione del 
target una volta fissato il numero delle ancore e la loro disposizione. Si sono anche mostrati i risultati 
di alcune simulazioni basate sui modelli esposti nel primo capitolo, eseguendo la localizzazione di un 
ipotetico target all’interno di una semplice costellazione di 4 ancore disposte ai vertici di un quadrato 
per fissare le idee. È lecito chiedersi se i risultati ottenuti dipendano dalla posizione del target rispetto 
alla posizione delle ancore, e in caso affermativo se esista un modo per quantificare a priori il livello 
di errore che ci si debba aspettare nel processo di stima della posizione. 

Quanto detto dipende dal fatto che ad ogni metodo affrontato può essere associata una differente 
interpretazione geometrica, caratterizzata da differenti famiglie di curve. A parità di perturbazione 
sull’errore iniziale, infatti, le intersezioni delle varie curve non sono caratterizzate dalla stessa pro- 
pagazione dell’errore andando ad agire in regioni diverse dello spazio disponibile; ciò è evidente se si 
pensa all’intersezione di rette piuttosto che a quella di cerchi. 

Nel capitolo precedente si sono testate le varie implementazioni dei metodi di localizzazione ope- 
rando sulla stessa posizione di riferimento, scelta in modo del tutto arbitrario. Si è osservato come 
tutti gli algoritmi scelti si comportino in modo simile operando su tale posizione, e si è apprezzato 
come tutti presentino lo stesso andamento asintotico nei confronti del rapporto segnale-rumore. Per 
valutare se questi risultati dipendano o meno dalla posizione del target si procede qui ad un’anali- 
si globale, ovvero si esegue l’algoritmo di localizzazione campionando l’intera area monitorata dalle 
ancore. Questo tipo di analisi si traduce nello spostare il target in maniera regolare lungo tutto il 
dominio. Operativamente si utilizza il passo costante Ax = Ay = 1m. 

In Figura 3.1 sono riportati i risultati di due esperimenti di triangolazione eseguiti con il metodo 
LLS introdotto nella sezione 2.1.1. In entrambi i casi si sono utilizzate 10 medie per ogni posizione e si è 
fissato il valore del rapporto segnali rumore a 30 dB. Per ogni posizione in cui si trova effettivamente il 
target è stata tracciata una freccia fino alla posizione stimata, pertanto a frecce più brevi corrispondono 
stime più accurate. 

I due quadri sinottici mostrati a destra e sinistra in Figura 3.1 sono formalmente identici come 
procedura di localizzazione e condizioni al contorno, ma del tutto scorrelati nei confronti del rumore: 
in entrambi i casi si può osservare come gli errori siano più contenuti in alcune zone (prevalentemente 
in un intorno del centro del dominio), mentre tendono ad aumentare spostandosi lungo il bordo del 
dominio. Dal momento che non c’è alcuna correlazione o effetto di memoria tra le due simulazioni, è 
ragionevole aspettarsi che la distribuzione dell’errore dipenda da qualche caratteristica intrinseca al 
metodo di localizzazione. 

Figura 3.2 descrive la stessa situazione in termini di dominio, posizioni del target e condizioni 
al contorno; in questo caso però l’algoritmo impiegato è il WLLS a due step, una delle migliori 
implementazioni possibili per la triangolazione, introdotto nella sezione 2.1.2. 

Benché la distribuzione dell’errore in questo caso sia diversa da quella descritta in Figura 3.1, si è 
mantenuto comunque il suo andamento nei confronti del dominio, ovvero gli errori minori si osservano 
verso il centro della costellazione. 

È evidente che se si portasse il numero di misure da 10 ad un valore asintoticamente alto, l’errore 
medio nella singola posizione si attenuerebbe fino a diventare nullo, dato che per tutti i metodi 
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Figura 3.1: Risultati di due esperimenti di localizzazione estesi ad un grigliato regolare di posizioni del 
target. In entrambi i casi l'algoritmo utilizzato & la triangolazione di tipo LLS. 
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Figura 3.2: Risultati di due esperimenti di localizzazione estesi ad un grigliato regolare di posizioni del 
target. In entrambi i casi l'algoritmo utilizzato è la triangolazione di tipo WLLS a due passi. 


presentati si è supposta l’assenza di errori sistematici. Questi ultimi infatti, qualora presenti sarebbero 
da imputare alla tecnica di misura del range — o dell’angolo di arrivo. Eventualmente sarebbe possibili 
eliminarli o al limite mitigarli con una procedura di calibrazione applicata al sistema di misura; tuttavia 
la questione va al di là dello scopo del presente testo. 

È invece interessante valutare lo scarto quadratico dell’errore in senso asintotico. Per stimarlo è 
necessario procedere a un numero elevato di test di localizzazione e valutare per ognuno di esso l’errore 
commesso nella stima delle coordinate x e y in valore assoluto. Il risultato di questa valutazione è 
presentato in Figura 3.3. Osservando la figura diventa evidente come la distribuzione dell’errore non 
sia casuale ma segua bensì un preciso andamento. Non è difficile concludere che i punti centrali 
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sono avvantaggiati perché equidistanti dalla costellazione delle ancore, per cui gli errori casuali che 
avvengono nella stima della posizione mediamente tendono ad elidersi, essendo il rumore a media 
nulla. 
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Figura 3.3: errore quadratico medio ottenuto stimando la posizione del target utilizzando l'algoritmo di 
triangolazione LLS introdotto in sezione 2.1.1. Le 91 posizioni uniformemente distribuite nel dominio sono 
state stimate 10000 volte per valutare l'errore quadratico medio. 


La stessa analisi eseguita sull’algoritmo WLLS a due passi è presentata in Figura 3.4. In questo 
caso la distribuzione dell’errore rispetta la stessa osservazione emersa nel caso precedente, ovvero una 
presenza di errori più bassi al centro. Tuttavia le due distribuzioni non sono esattamente sovrapponibili: 
si nota infatti che, a differenza di quanto avviene in Figura 3.3, in questo caso l’errore si mantiene 
relativamente contenuto anche sulla congiungente di due ancore, lungo i lati del quadrato. Si può 
quindi fin qui concludere che una causa della distribuzione dell’errore sia dovuta alla disposizione 
delle ancore (data la similitudine dei risultati ottenuti con due metodi diversi), ma che al contempo 
esista comunque una componente di questa distribuzione indubbiamente figlia del metodo scelto. 


3.1 Il limite teorico della precisione: il Cramer Rao Bound 


Lo strumento analitico utilizzato per affrontare questo tipo di analisi è la teoria della quantità di 
informazione sviluppata da Fisher, e in particolare uno dei suoi risultati più notevoli, il criterio noto 
come Cramer Rao Bound, CRB, talvolta indicato anche come Cramer Rao Lower Bound, CRLB, [36]. 
Questo indice metrico fornisce una stima della minor varianza teoricamente ottenibile una volta fissato 
uno stimatore, nell’ipotesi che sia nota la statistica di rumore dei dati in ingresso. Per motivi pratici 
in genere ci si riduce a considerare soltanto statistiche gaussiane dell’errore, data la difficoltà insita nel 
modellizzare densità di probabilità di tipo qualsiasi. Benché questa ipotesi possa sembrare restrittiva, 
data la grande quantità di processi e meccanismi che a seconda del metodo possono colorare la statistica 
dell’errore, l’approccio qui descritto trova comunque grande riscontro anche nella verifica sperimentale. 
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Figura 3.4: errore quadratico medio ottenuto stimando la posizione del target utilizzando l'algoritmo di 
triangolazione LLS introdotto in sezione 2.1.2. Sono state considerate le stesse posizioni con le stesse 
condizioni di Figura 3.3. 


Il CRB si configura come l’asintoto per le prestazioni dei vari metodi di localizzazione, ed è in 
generale una funzione del solo modello di osservazione, non dello specifico metodo risolutivo. Si può 
anzi utilizzare come metro di paragone di ogni algoritmo, dato che l’efficienza di un approccio può 
essere misurata valutandola in rapporto a quanto previsto dal CRB nella medesima configurazione. 
In quanto varianza limite, il valore di CRB va confrontato con lo scarto quadratico medio dell’errore 
di posizionamento, e pertanto fornisce una valutazione aprioristica della precisione del metodo in 
esame. È bene precisare infatti che tutti i metodi qui valutati sono non polarizzati (unbiased), ovvero 
forniscono in media il valore esatto. Se esiste un bias, infatti, esso si traduce in un errore sistematico 
delle misure, e può pertanto essere sempre corretto. 

Il punto di partenza per la costruzione del CRB è la matrice dell’informazione di Fisher, nota con 
il nome di Fisher Information Matrix (FIM). Quest'ultima è una metrica che valuta la quantità di 
informazione contenuta in un modello di osservazione. 

I passi per calcolare il CRB si possono schematizzare nella seguente procedura: 


1. si calcolano le derivate seconde del logaritmo della distribuzione di probabilità della misura 
0° In p(r). 
(929%) 


rispetto all’incognita posizionale, ovvero: 


2. a partire da queste grandezze, si costruisce la matrice di Fisher stimando il valore atteso dei 
valori calcolati al punto 1, ovvero: 


o? a (3.1) 


PAS er 


con F(x) che indica la matrice di Fisher; 


3. DILUIZIONE DELLA PRECISIONE 45 


3. il valore del CRB si ottiene a questo punto invertendo F (æ). In particolare, nel caso in questione, 
la funzione CRB per gli elementi in x e y sono le componenti diagonali dell'inverso della matrice 
di Fisher, ovvero: 


II 


CRB(x) =[F (x)] 
CRB(y) =[F (æ)] 


(3.2a) 
(3.2b) 


SLI] 
NE ee 


Come anticipato, nel caso in cui le misure siano affette da un errore gaussiano a media nulla, la 
matrice di Fisher può essere calcolata in modo ancor più agevole, nel seguente modo: 


n= [eT e [4 


dove C' è la matrice delle covarianze per il rumore n. A partire da quest’ultima espressione è possibile 
valutare in maniera diretta il CRB dei vari metodi di localizzazione. Si ricorda ancora come questo 
dipenda dal modello di osservazione, comprendendo in questo la specifica morfologia della costellazione 
di ancore, ovvero disposizione e densità nei confronti del dominio monitorato. 


(3.3) 


3.2 Cramer Rao Bound per la trilaterazione 


Si procede quindi a stimare la forma della funzione CRB per l’approccio di localizzazione basato 
sulla trilaterazione. Ripartendo dalla funzione di osservazione in questione: 


V (a= 21)? + (y=) 
V (2 = 22)? + (y — y2)? 

f(x) = (3.4) 
v (zu)? + (y — yr)? 


? 


è semplice realizzare che il primo passo per il calcolo del CRB, ovvero la derivata direzionale di (3.4), 
risulta: 


L_X1 Y Y1 
V(&-21)2+(y-y1)? V(&-21)2+(y-y1)? = ra 
v2 y—y2 2-22 y—y2 
Of p(x) _ | Vla-z2)>+(y-y2)? y (e—=2)2+(y=ya)? | _ da d2 35 
Ox u : . = : : ’ ( y ) 
z-s, v-yL TIL Y YL 
A Vee +U) i a 
dove si è sfruttato la definizione di distanza tra target e costellazione: 
V (w= 21)? + (y =y)? = d. (3.6) 


Ricordando la definizione della matrice di correlazione e sostituendo l’espressione (3.5) in quella 
semplificata riportata in (3.3) si ottiene: 


2 
o 0 0 2-21 YY 
R,1 dí di 
e al È L-L2 y Yy2 
T—L1 TX L_XL 0 
o d d 
= di da dí R,1 2 2 
Fr(&) = as, y—y2 a) > (3.7) 
dy da dí 
X_XL Y-YL 
0 ORL de di 
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calcolando la quale si perviene alla forma esplicita della matrice di Fisher: 


L 5 L 
a ar 
Fra)=| ¿2% FR | (3.8) 
(2-21) (y-41) (v=yı) 
2 orid 2 RA 


CRB,(2) = [Fra], = LI (3.9a) 


—1 l > 
CRB, (a) = [Pr(e)), > = >> (3.9b) 
con |A| determinante della matrice A. 
Da queste, si puó pervenire alla varianza assoluta sulla localizzazione, ovvero alla varianza dell'er- 


rore di localizzazione del tipo |£ — x|, come somma delle due componenti 
CRBr(x)=CRB,(x)+CRB,(x). (3.10) 


Come anticipato, questa funzione non contiene alcuna informazione sullo specifico metodo risolutivo, 
ma soltanto l’effetto della costellazione sull’area monitorata, dal momento che ogni punto dello spazio 
osservato dal CRB viene messo in relazione all’insieme delle distanze relative dalle ancore attraverso le 
distanze d;. Partendo dagli altri metodi di localizzazione, basati su modelli di osservazione diversi, si 
ottengono espressioni simili ma non identiche. Nel caso della triangolazione differenziale, per esempio, 


; ə . 
la matrice Ea risulta: 
T—T2 L_X1 y y2 y Yi 
da 1 2 1 
L—L3 221 Y—Y3 Y=zyı 
d3 di d3 dı 
For(®) = A E ; (3.11) 
LLL L—L1 Y=YL y=yı 
di di dL di 
per la trilaterazione basata su RSSI: 
x— Tı Y=yı 
dí dí 
T-T2 Y=y2 
a2 d2 


IRrssı(®) = -a : ; (3.12) 


per la triangolazione: 


y—yı L_X1 
di dî 
— Ya 
di 
Fılz) = E bh (3.13) 
your tn 
di, di 


e così via per ogni modello di osservazione. 
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3.2.1 Alcuni esempi di CRB per costellazioni semplici 


Per valutare l’effetto della costellazione sulle prestazioni che ci si può attendere da un qualsiasi 
algoritmo di localizzazione si illustra qui l’effetto del CRB in accordo con le equazioni (3.9). Con 
riferimento alla Figura 3.5, il grafico a sinistra implementa il CRB per la coordinata x secondo l’equa- 
zione (3.9a), mentre il grafico a destra utilizza l’equazione (3.9b) per la coordinata y. Un’informazione 
riassuntiva si può avere osservando invece la Figura 3.6 che mostra i risultati del CRB totale come da 
equazione (3.10). 


10 0.9 10 x 5 0.9 

E 0.8 E A 0.8 

0.7 Koa 0.7 
E? 0.64 Et 0.6 © 
> 5 > Pi 
jo} O ES 
a 4 0.55 34 0.55 

0.4 — > 0.4 

È 0.3 £ NY 0.3 

ô 0.2 0 x o } 0.2 

1 0.1 

0 2 4 6 8 10 9 0 2 4 6 8 10 
lato x (m) lato x (m) 


Figura 3.5: CRB per la coordinata x, a sinistra, e y, a destra, in accordo con la formula (3.9). 


Confrontando la Figura 3.6 con la Figura 3.3 e la Figura 3.4 si conferma quanto anticipato nel 
paragrafo precedente: l'errore quadratico medio dell'algoritmo LLS in Figura 3.3 e WLLS a due passi 
in Figura 3.4 ricalca la struttura prevista dalla figura Figura 3.6. In particolare, si ha la conferma della 
qualitá del metodo WLLS a due passi, il quale segue in maniera molto precisa il profilo previsto dalla 
funzione CRB (come si può osservare confrontando le regioni a metà dei lati del quadrato che definisce 
la costellazione). Si può affermare che questo algoritmo è asintoticamente efficiente, dal momento che 
ottiene le prestazioni limite previste dalla teoria dell’informazione. Si può quindi definire un’efficienza 
di localizzazione generalizzando il concetto con la seguente forma: 


E [(@ — x)? + G- yy] 
CRB(x, y) i 


n(z, y) = (3.14) 
dove le stime F e Y sono ottenute con lo specifico algoritmo di cui si sta valutando l’efficienza. 

Dal momento che non si può ottenere un risultato migliore di quello previsto dal CRB (almeno 
nel caso di algoritmi non polarizzati), se le prestazioni ottenute non sono adeguate agli obiettivi dello 
specifico problema di localizzazione, l’unica alternativa è modificare la costellazione delle ancore. È 
infatti facile comprendere come il numero e la disposizione delle ancore sia infatti il punto chiave per 
ottenere prestazioni asintotiche migliori, dal momento che il denominatore delle (3.9) cresce con L più 
rapidamente di quanto faccia il numeratore. 

In Figura 3.8 è mostrato il risultato ottenuto aggiungendo una singola ancora al centro del dominio 
da monitorare. La distribuzione della varianza limite si è chiaramente modificata di conseguenza e in 
particolare, confrontando Figura 3.8 con la corrispondente in Figura 3.6 si nota come il valore minimo 
sia molto diminuito. Elevati valori di CRB sono comunque presenti nelle regioni adiacenti alle boe: 
questo succede perché la propagazione di errore continua ad essere sfavorevole nonostante la densità 
di ancore. 
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lato x (m) 


Figura 3.6: CRB complessivo per la costellazione di 4 ancore, in accordo con (3.10). Questa distribuzione & 
valida per qualsiasi approccio basato sulla triangolazione che risponda al modello di osservazione in (1.12), 


indipendentemente dalla specifica implementazione utilizzata (lineare, pesata, non lineare, etc.). 


IKON 


10 


lato y (m) 


Figura 3.7: CRB per la coordinata x, a sinistra, e y, a destra, in accordo con (3.9), quando alla costellazione 
iniziale descritta in Figura 3.5 viene aggiunta una quinta ancora al centro del dominio. 


Una costellazione di 8 elementi produce i risultati mostrati in Figura 3.9 e in in Figura 3.10. 
Con questa configurazione si ottengono prestazioni asintotiche senza dubbio migliori, ma il prezzo da 
pagare è una densità di ancore per superficie monitorata molto elevata. 

Modificando la posizione degli elementi della costellazione si può ottenere un risultato migliore 
senza ricorrere a eccessive densità delle ancore. Con riferimento alla Figura 3.11, 4 boe sono state 
posizionate sui punti medi dei lati piuttosto che sui vertici del quadrato. Si può osservare che, nell’area 
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lato x (m) 


Figura 3.8: CRB complessivo per la costellazione di 5 ancore, in accordo con la (3.10). 
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Figura 3.9: CRB per la coordinata x, a sinistra, e y, a destra, in accordo con (3.9), operando su una 
costellazione di 8 ancore distribuite lungo il contorno dell’area monitorata. 


sottesa dalle 4 ancore, la distribuzione dell’errore è identica a quella del caso Figura 3.6, ma potendo 
contare su distanze ridotte di V2, la distribuzione dell’errore raggiunge al massimo la metà di quello del 
primo caso. Inoltre, anche nelle zone periferiche dell’area osservata l’errore non cresce eccessivamente. 
In definitiva la situazione risulta molto vantaggiosa se misurata con la metrica del massimo CRB per 
numero di ancore e per area monitorata. 
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lato x (m) 


Figura 3.10: CRB complessivo per la costellazione di 8 ancore, in accordo con (3.10). 


lato x (m) 


Figura 3.11: CRB di una costellazione di 4 ancore posizione sui punti medi del bordo dell'area monitorata. 


3.3 Analisi di media e varianza: l’approccio ottimo 


Nei casi analizzati, ogni metodo di stima posizionale si riconduce alla massimizzazione (o mini- 
mizzazione) di una funzione continua. In questo contesto, & possibile valutare la media e la varianza 
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della posizione stimata © utilizzando gli strumenti dello sviluppo in serie nell'intorno del massimo (o 
minimo) del relativo funzionale. Si indichi con J(#) una generica funzione continua di € e ® sia il suo 
massimo (o minimo), [37]. In queste condizioni si può certamente affermare che: 


aJ(#) 
di 


= 0. (3.15) 


Ipotizzando che, in presenza di errori di stima contenuti, £ sia ragionevolmente vicino alla soluzione 
esatta x, l’espansione di Taylor attorno a questo valore fornisce la seguente espressione: 


OJ (x) OJ (2) 0? J(£) = 

0 = = ® — H —2). 3.16 
ee ern ree 

Da questa uguaglianza, ancorché approssimata, si desume che 

OJ (2) 0? J(T) Sa 
_ — x TI), 3.17 
da |: 2508 lex ( ) (3.17) 
che può essere convenientemente riscritta come 

-V(J(a)) = H(J(e))(£- x), (3.18) 


dove H(J(x)) e V(J(x)) sono rispettivamente la matrice Hessiana e il vettore gradiente del funzionale 
J(x), valutati nella posizione reale. Nel caso in cui le derivate seconde della matrice Hessiana siano 
abbastanza regolari nell’intorno di x, si può porre, [38]: 


H(J(a)) = E[H(J(2))]. (3.19) 


Utilizzando la (3.15) e valutando il valore medio dei due membri di (3.16), si perviene alla stima del 
valore medio di X, che assume la seguente forma: 


E [2] =x-[EH(J())] t E[V (J (æ))]. (3.20) 
Se Z è una stima non polarizzata — ovvero unbiased — di x, si ha che E [£] = a, e di conseguenza 
Pultimo termine in (3.20) è pari a zero. Utilizzando di nuovo le equazioni (3.18) e (3.19) e le proprietà 
di simmetria delle matrice Hessiana, si può determinare la covarianza di €, qui indicata con il nome 
di Cz: 
C.=E[@-El@)@-E@)]=E[@-o)@-2)]= 
[E [H (J (æ) E [VJ (æ)V" J(æ)] EH). (3.21) 


A partire da questa matrice, le varianze delle stime di x e y sono date dalle componenti diagonali 
[Celli e [Calo 2: 


3.3.1 Analisi del metodo Maximum Likelihood 


Alla luce dell’analisi del paragrafo precedente, si riprende adesso in considerazione il metodo 
di localizzazione ML (massima verosimiglianza), espresso dall'equazione (2.96). Dal momento che 
E [rr] = di, il valore atteso della matrice Hessiana può essere facilmente calcolato a partire dalle 
equazioni (2.100), e risulta pari a: 


CRI Tr, 


L L 

(2-21) (y—Y1) (y=yi)? 
5 o? d2 y 2 d2 
l=1 E11 l=1 


L LL: 2 L = = 
2 ee x E 


E[H(Jur,r(2))] © 2 (3.22) 


RIU 
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inoltre, il valore atteso del gradiente risulta: 
E[V(Imz,r(x))] = 0. (3.23) 


Dal momento che H(Jur,r(&)) è non singolare, è chiaro dall’equazione (3.20) che la stima di tipo ML 
è non polarizzata. Inoltre, considerando l’equazione (2.101) e ricordando che, per definizione, risulta 


E Fr = ay] = Rn si ha: 


Se Y e) 
ML,R(® x)| = = = . 
E [Jmr R(®) Jit R i el a ET 3.24 
> (2-20) (Y-41) y (uu)? 
í=1 rd í=1 Rid 
Sostituendo le uguaglianze in (3.23) e (3.24) nell’equazione (3.21) si osserva che: 
Y gw y E) 
A | p5 At ART |= TRæ) (3.25) 
$ (2-20) (yy) EN (y=)? 
í=1 Toa LY í=1 Foa 1t 


Quest'ultima uguaglianza dimostra come lo stimatore di tipo Maximum Likelihood sia ottimo, nel 
senso che la sua varianza coincide con il valore del CRB. In altre parole l’approccio ML satura tutta 
l'informazione disponibile ottenendo la migliore stima possibile in senso statistico, nel caso in cui il 
rumore segua una distribuzione gaussiana a valor medio nullo. 


3.3.2 Analisi del metodo Weighted Linear Least Square 


La stessa analisi di media e varianza può essere applicata anche agli approcci lineari. Considerando 
la funzione di costo di tipo pesato, ovvero il funzionale WLLS definito dall’equazione (2.53) e il 
suo successivo svolgimento, si procede anche in questo caso con il definire la corrispondente matrice 
hessiana: 


OJwLLs,r(P) T 
A ARA =H = 2A A 2 
000pT _ (Jwrrs,r(P)) W > (3 6) 
e quindi il vettore gradiente: 
Ms smtp) = 2A" W (Ap —b). (3.27) 
P=p 


Sfruttando i risultati delle equazioni e (2.12) e (2.17), il valore atteso del gradiente V(Jwrrs,r(0)) 
risulta: 

E[V (Jwrrs.r(p))] = 2A' W E [q] = 0. (3.28) 
Dalle equazioni (3.20), (3.26) e (3.28) si può evincere come il metodo WLLS fornisca stime non 
polarizzate per la posizione x. Inoltre, utilizzando l’equazione (2.56) è possibile calcolare il valore 
dell’espressione V(Jwrrs,.r(pP))V (JwLLs,R(p)) in forma esplicita, che risulta: 


E [V Jwrrs,r(P)) Y" (Jwrrs,r(p))] = ¿A WE faq] WA = ¿A WA; (3.29) 


sostituendo l’equazione (3.26) e l’equazione (3.28) nell’equazione (3.21), si ottiene la matrice di 
covarianza per P, indicata con Cp, pari a: 


C, = (A'wA) — l (3.30) 


Questa espressione è in perfetto accordo con l’analisi prevista dalla teoria BLUE [32] nella forma 
dell’equazione (2.65). 
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3.4 Considerazioni finali sulla costellazione di ancore 


Per riassumere i vari aspetti discussi in questo capitolo, in figura Figura 3.12 è descritto il setup 
di un esperimento di localizzazione simulato caratterizzato da una costellazione variabile di ancore. 
Utilizzando questo schema, si vogliono mettere a confronto le prestazioni dei vari algoritmi con il limite 
sancito dal CRB al variare della configurazione e della densità delle ancore. In paticolare, fissato il 
punto æ = (4, 1), si eseguono i vari algoritmi di localizzazione aumentando progressivamente le ancore 
attive nella costellazione. 
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Figura 3.12: Costellazione variabile delle ancore. 


Fissato il rapporto segnale rumore al valore di 30 dB, la Figura 3.13 mostra i risultati ottenuti con 
gli algoritmi LLS, WLLS, WLLS a 2 fasi, NLS di tipo Newton-Raphson e ML di tipo Newton-Raphson. 
Le varie stime sono ottenute passando da L = 3, il minimo per identificare la posizione, fino a L = 25. 

Per ogni costellazione è stato calcolato il corrispondente valore del CRB. Osservando la figura si 
può apprezzare come, ancora una volta, sono gli algoritmi Maximum Likelyhood e Weighted Least 
Square a due fasi a eguagliare il limite asintotico del CRB. Si conferma altresì che l'approccio LLS 
ha le peggiori prestazioni a parità di condizioni iniziali. È inoltre interessante come algoritmi diversi 
reagiscano diversamente alla disposizione delle ancore. Per esempio, l'algoritmo WLLS ha un netto 
miglioramento all'introduzione della 17”°® e 18" ancora, mentre l’algoritmo NLS è poco influenzato 
da questo incremento di densità. 

La stessa analisi è stata compiuta di nuovo, andando a posizionare però casualmente le 25 ancore 
nella stessa area. La disposizione risultante è mostrata in Figura 3.14, mentre i relativi risultati in 
Figura 3.15. Data la natura casuale delle posizioni delle ancore è più difficile trarre conseguenze 
sull’ottima costellazione delle stesse, ma allo stesso tempo è ancor più neutrale il giudizio sulla qualità 
dei metodi risolutivi. Analizzando quindi i risultati si conferma ancora una volta come il metodo ML 
e il metodo WLLS a due passi producano risultati paragonabili al limite del CRB. 
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Figura 3.14: Costellazione variabile delle ancore. 
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Figura 3.15: Costellazione variabile delle ancore. 
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GUIDA AL MATERIALE DI APPROFONDIMENTO 


In questo capitolo viene fornita una guida al materiale di autoapprendimento sugli innumerevoli 
temi implicati nel mondo della localizzazione. Questo elenco non ha la pretesa di essere esaustivo, e si 
propone come riferimento per il lettore in modo che possa approfondire i moltissimi aspetti solamente 
accennati nel presente testo, che, come dichiarato fin dall’inizio, aveva come scopo quello di essere 
soltanto una guida introduttiva a un mondo estremamente vasto e in costante divenire. 


— S. Sand, A. Dammann, C. Mensing, Positioning in Wireless Communications Systems. Wiley & 
Sons, 2014. 


Questo testo e incentrato sulle principali differenze e somiglianze tra sistemi di comunicazione wireless 
e sistemi di navigazione. Vengono discussi scenari che sono critici per i sistemi di navigazione, come il 
sistema di posizionamento globale (GPS), e che motivano l'uso di sistemi di posizionamento basati sui 
sistemi di comunicazione wireless terrestri. Il libro introduce diversi approcci per la determinazione dei 
parametri che dipendono dalla posizione del terminale mobile e offre una disamina di algoritmi per la 
stima della posizione del target. Vengono inoltre proposti modelli per il canale radio e la mobilitá degli 
utenti, aspetti importanti da tenere sotto controllo nel caso di esperimenti in ambienti reali e dunque 
per implementare simulazioni piú indicative. Scritto da esperti del settore che lavorano all'avanguardia 
dello stesso, gli autori forniscono una visione eccellente su questo argomento, consentendo una copertura 
approfondita degli attuali sviluppi in questo campo in rapida espansione. 


— R. Mannings, Ubiquitous Positioning. Artech House, 2008. 


Questa testo innovativo è un’ottima introduzione per il campo emergente della consapevolezza spaziale, 
con speciale riferimento al campo sensoristico, dove sistemi di posizionamento che permettono di loca- 
lizzare oggetti, veicoli o persone stanno diventando realtà. Il futuro e la crescita di questa tecnologia 
risulterà dalla convergenza di molti altri settori, tra cui la telematica mobile, la tecnologia dell’ Internet 
of Things, il web semantico e così via. Questo testo, prima del suo genere, il volume lungimirante, che 
esplora il mondo della consapevolezza spaziale da una prospettiva di convergenza, offrendo una road 
map per questo settore in rapida crescita. 


— T. E. Tuncer, B. Friedlander, Classical and Modern Direction-of-Arrival Estimation. Academic 
Press, 2009. 


Il testo, scritto da esperti del settore, presenta in maniera esaustiva la teoria e la pratica della stima della 
direzione di arrivo di segnali arbitrari, fondamento delle tecniche di triangolazione. Questa commistione 
unica di tecniche utilizzate nell’ambito di sistemi commerciali e di metodi super-risolutivi allo stato 
dell’arte è una preziosissima fonte di informazioni per ingegneri e ricercatori che si vogliono cimentare 
nel campo della localizzazione indoor. 


I principali argomenti trattati sono: metodi classici per la stima della direzione di arrivo, soluzioni 
adottate nei sistemi commerciali, calibrazione di sistemi di array di antenne, algoritmi per segnali a 
banda larga, distribuzione della direzione d’arrivo nel dominio tempo-frequenza, stima della direzione 
di arrivo a soglia, localizzazione con reti di sensori. 


— Alan Bensky, Wireless positioning technologies and applications, Artech House, 2016. 
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Nonostante il successo del GPS, sono i metodi basati su comunicazioni a corto raggio ad avere davanti 
a se la maggiore prospettiva di crescita, e in particolare le reti di sensori per applicazioni industriali e 
domotiche. 


Questo testo raccoglie tutti questi metodi insieme fornendo un riferimento unificato, analizzando van- 
taggi e svantaggi, e le procedure di attuazione per tutti i principi di funzionamento considerati. Si tratta 
di un valido aiuto per la scelta della tecnologia per chi si affaccia per la prima volta in questo mondo, 
aiutando il lettore a soppesare i costi in relazione ai migliori risultati possibili. Questa risorsa copre i 
principi fondamentali dei metodi basati sul tempo di volo (ToA, TDoA), i metodo basati sulla potenza 
del segnale (RSSI), e quelli basati sull’angolo di arrivo (DoA). 


S.A. Zekavat, R. M. Buehrer, Handbook of Position Location, Theory, Practice, and Advances. 
Wiley, 2011. 


I sistemi radio in grado di offrire una stima della posizione dell’utente rappresentano un’applicazione 
emergente in molti campi, tra i quali sicurezza e monitoraggio ambientale, sorveglianza e applicazione 
della legge, interventi di emergenza e di protezione civile, coordinamento interveicolare per il rilevamento 
del rischio di collisione e in molti altri campi ancora. 


Sebbene la rapida ascesa di queste tematiche, gli ambienti universitari e accademici sono tuttora impre- 
parati alla divulgazione delle conoscenze teoriche e degli strumenti necessari alla base di queste moderne 
applicazioni. Il testo in esame si propone come guida per colmare questo vuoto tra le richieste sempre 
più attuali del mercato e la carenza di informazioni relativi a questo campo. In particolare, il testo si 
propone come manuale adatto all'autoapprendimento, fornendo un ampia scelta di algorirmi e modelli 
implementati in MATLAB per una rapida fruizione degli argomenti trattati. 


S. Chandran, Advances in Direction-of-Arrival Estimation. Artech House, 2005. 


Determinare la direzione di arrivo (DOA) di un segnale arbitrario è di importanza critica per gli inge- 
gneri elettronici ed elettromagnetici che lavorano nel campo delle applicazioni civili ma anche nel campo 
militare. L'importanza della direzione di arrivo é insita nel fatto che su questa grandezza si fondano i si- 
stemi di triangolazione, i quali stanno rapidamente soppiantando i sistemi di trilaterazione, specialmente 
in ambito radar per la rilevazione rapida di bersagli. Questa testo all’avanguardia presenta gli ultimi 
sviluppi nella stima dell'angolo di arrivo presentati dai principali ricercatori operanti nel campo. Il testo 
e un riferimento approfondito per i moderni approcci operativi e copre una vasta gamma di argomenti 
chiave, da metodi di stima alle applicazioni specifiche di stima DOA e la configurazione sperimentale e 
risultati. 


James Bao-Yen Tsui, Fundamentals of global positioning system receivers: a software approach. 
John Wiley & Sons, 2005. 


Sebbene complementare ai temi trattati nel presente testo, la conoscenza del sistema GPS risulta co- 
munque fondamentale per comprendere la necessitá degli approcci indoor. Questo testo consente ai 
lettori di comprendere e applicare i complessi principi di funzionamento del sistema di posizionamento 
GPS. L’autore, una delle principali autoritá nel campo GPS, presenta il materiale secondo l’approccio 
della software defined radio applicato al ricevitore GPS, un approccio che aiuta i lettori a comprendere 
meglio il funzionamento e che riflette l'integrazione del ricevitore nei dispositivi di uso quotidiano, come 
i telefoni cellulari. 


Il libro fornisce gli strumenti e le informazioni necessarie per comprendere e sfruttare tutti gli aspetti 
della tecnologia dei ricevitori, nonché schemi di navigazione rilevanti: panoramica delle basi GPS e la 
costellazione di satelliti che compongono il sistema GPS, esame dettagliato dei segnali coinvolti e loro 
acquisizione e contestuale tracciamento, presentazione e analisi degli algoritmi necessari a calcolare la 
posizione sulla base dei segnali acquisiti, dimostrativo software del sistema dei passi descritti. 


— Steven M. Kay, Fundamentals of Statistical Signal Processing, Prentice Hall, 2013. 
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In questo testo, l’autore, uno dei massimi esperti mondiali sull’argomento, offre una delle pit complete 
panoramiche sul rilevamento e al stima del segnale, sulla sua elaborazione e sullo sviluppo di algoritmi 
efficienti, coprendo un ampio spettro di aspetti, dalla propagazione dell’errore e relativa simulazioni 
al computer, fino alla valutazione delle prestazioni. L’autore propone strategie e algoritmi che hanno 
“superato la prova del tempo,” e offre esempi reali mutuati da diverse aree applicative. 


Gli argomenti trattati in questo testo, che si compone di più volumi, comprendono: approccio graduale 
alla progettazione di algoritmi, il confronto e la scelta di modelli di segnale e rumore, la valutazione delle 
prestazioni, gli approcci ottimali utilizzando i “grandi teoremi”, gli algoritmi per la stima, la rilevazione 
e stima spettrale. Infine l’autore propone casi di studio completi: il radar doppler, il rilevamento del 
segnale magnetico, e il monitoraggio della frequenza cardiaca. 
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Introduzione ai Sistemi di Localizzazione 
Indoor costituisce il secondo numero della 
collana «Quaderni del Dottorato di Ricerca in 
Ingegneria dell'Informazione» dell'Università 
di Firenze, a supporto didattico dell'omonimo 
corso di dottorato. Questo quaderno è pensato 
come testo introduttivo per chi voglia cimentarsi 
nel campo emergente della localizzazione di 
mezzi mobili in ambienti indoor, dove i servizi 
di posizionamento satellitari sono assenti o 
inaffidabili. Nello spazio di quattro capitoli, 

il presente testo passa in rassegna i metodi 
principali di localizzazione, i corrispondenti 
algoritmi, e la propagazione dell'errore 
associata a questi ultimi. 
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alle antenne intelligenti per applicazioni di 
localizzazione indoor. 


